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ABSTRACT 

An  analytical  and  numerical  analysis  was  carried  out  on  the  operation  of 
an  electrodynamic  tether  system  used  for  power  generation  and/or  thrusting 
in  a  space  environment.  Three  problems  were  examined.  First,  the 
efficiency  of  an  uninsulated  tether  of  prescribed  design,  as  determined  by 
the  magnitude  of  current  leakage  along  its  length  due  to  positive  ion 
capture  and  secondary  electron  emissions,  was  compared  to  that  of  a 
perfectly  insulated  tether  of  identical  design.  Second,  the  effects  on 
system  efficiency  of  variations  in  the  design  parameters  of  the  tether  and 
the  orbit  in  which  it  operates  were  examined  by  means  of  a  numerical 
analysis  of  the  thermal  balance  of  the  system.  Third,  the  effects  which 
the  mode  of  operation  of  the  tether  has  on  the  classical  elements  of  the 
orbit  in  which  it  operates  were  examined  through  a  numerical  analysis. 

The  results  of  the  first  examination  indicated  that  the  uninsulated  tether 
can  still  operate  as  a  power  generator  or  thruster,  although  with  an 
effective  efficiency  when  compared  to  the  perfectly  insulated  tether  of 
about  90%.  The  second  examination  indicated  that  tethers  operated  more 
efficiently  at  lower  temperatures,  and  that  these  temperatures  were 
directly  affected  by  the  surface  characteristics  of  the  tether, 
specifically  the  absorptance  and  emittance.  High  emittance  and  low 
absorptance  produced  the  lowest  average  operating  temperatures,  and  hence 
the  highest  efficiencies.  Fluctuations  in  the  temperatures  over  an  orbit 
were  about  ±  20  degrees.  These  variations  were  generally  unaffected  by 
the  orientation  of  the  tether's  orbit,  and  by  the  position  of  the  earth, 
with  respect  to  the  sun.  The  third  examination  indicated  that  the 
classical  elements  of  the  tether's  orbit,  with  the  desired  exception  of 
the  semi -major  axis  which  represents  orbital  energy,  were  generally 
unaffected  when  the  tether  was  operated  in  a  single  mode,  but  dramatically 
increased  the  eccentricity  of  the  orbit  when  operated  in  a  mixed  mode 
within  an  orbit,  so  much  so  that  operation  in  this  matter  is  prohibitive. 


Thesis  Supervisor:  Dr.  Manuel  Martinez-Sanchez 

Title:  Associate  Professor  of  Aeronautics  &  Astronautics 


TABLE  OF  CONTENTS 


SECTION 


PAGE 


Abstract  .  1 

Table  of  Contents  .  2 

List  of  Figures  .  5 

1.  INTRODUCTION  .  6 

2.  THEORETICAL  BACKGROUND  .  8 

2.1  INDUCED  VOLTAGE  .  8 

2.2  EQUIVALENT  CLOSURE  CIRCUIT  .  10 

2.2.1  Ionospheric  Impedance  .  11 

2.2.2  Contactors  .  12 

2.2.3  Insulation  .  16 

2.2.4  Reduced  Circuit  .  18 

2.3  MODEL  FOR  INDUCED  VOLTAGE  .  18 

3.  METHODS  OF  ANALYSIS  .  22 

3.1  INSULATION  .  23 

3.2  METHOD  OF  THERMAL  ANALYSIS  .  30 

3.2.1  Thermal  Balance  .  31 

a.  Time  Constant 

b.  Assumption  of  Uniform  Tether 

3.2.2  Calculation  of  Components  of  Thermal  Balance .  33 


a.  Radiation 

b.  External  Flux 

1.  Solar 

2.  Albedo 

3.  Earth  Radiation 

4.  Aerodynamic 

c.  Conductor  Material 


3.2.3  Conductor  Material 


38 


TABLE  OF  CONTENTS  Continued 


SECTION 


PAGE 


3.3  ORBITAL  VARIATIONS  .  39 

3.3.1  Disturbing  Acceleration  Vector  .  39 


a.  Induced  Electromotive  Force 

b.  Aerodynamic  Force 

c.  Solar  Radiation 

d.  Tether  Libration 

e.  Earth  Oblateness  Effects 


3.3.2  Variation  Equations  .  45 

4.  APPLICATION  OF  METHODS  &  INTERPRETATION  OF  RESULTS  .  46 

4.1  THERMAL  EFFECTS  .  46 

4.1.1  Effects  of  Tether  Parameters  .  47 


a.  Absorption  &  Emittance 

b.  Conductor  Material 

c.  Diameter 

d.  Length 

4.1.2  Effects  of  Orbital  Parameters  .  51 

a.  Altitude 

b.  Inclination 

c.  Line  of  Nodes 

d.  Seasonal  Variations 


4.2  VARIATIONS  OF  ORBITAL  ELEMENTS  .  54 

4.2.1  Open  Circuit  Operation  .  55 

4.2.2  Pure  Mode  Operation  .  57 

4.2.3  Mixed  Mode  Operation  .  59 

5.  CONCLUSION  .  62 

5.1  SUMMARY  .  62 

5.2  RECOMMENDATIONS  .  63 


3 


TABLE  OF  CONTENTS  Continued 


SECTION 


FIGURES  .. 
APPENDICES 
REFERENCES 


LIST  OF  FIGURES 


1  Equivalent  Circuit  Diagram 

2  Contactor  Configurations 

3  Ion  Density  Fluctuations 

4  Simplified  Circuit  Diagram 

5A  Voltage  Variations:  Tilted  Dipole  Model 

5B  Voltage  Variations:  Octopole  Model 

6  Voltage  Variations:  Expanded  Model 

7  Momentum  Balance  for  Ion  Capture 

8  Current  Direction  in  Generator  Mode 

9A  System  Potential  Drops  -  Generator  Mode 
9B  System  Potential  Drops  -  Thruster  Mode 

10  Rotation  Angles 

11  Variation  in  Tether  Properties 

12  Conductance/Mass  Comparison 

13  e,  p,  h  Coordinate  System 

14  Eccentricity  Vector 

15  Tether  Property  Variations  Over  an  Orbit 
16A  Effects  of  Absorptance  &  Emittance 

16B  Typical  Dielectric  Properties 

17  Average  Operating  Temperature  for  Different  Conductors 

18A  Effects  of  Diameter  with  Constant  Current 
18B  Effects  of  Diameter  with  Constant  Load  Power 

19  Effects  of  Length  for  Constant  Load  Power 

20  Effects  of  Orbital  Altitude 

21  Temperature  Variations  due  to  Line  of  Nodes 

22  Temperature  Variations  due  to  Sun  Angle 

23  Open  Circuit  Operation  -  Orbital  Element  Variations 

24  Generator  Mode  Operation  -  Orbital  Element  Variations 

25  Thruster  Mode  Operation  -  Orbital  Element  Variations 

26  Mixed  Mode  Operation  -  Orbital  Element  Variations 

27  Eccentricity  with  Increasing  Power-to-Mass  Ratios 


1 .  INTRODUCTION 


Electrodynamic  tethers  are  essentially  long  conducting  wires  which, 
when  deployed  from  an  orbiting  satellite,  can  generate  power  by  converting 
orbital  mechanical  energy  into  electrical  energy  for  use  on  the  satellite. 
They  can  accomplish  this  as  a  consequence  of  the  electrodynamically 
induced  voltage  caused  by  the  movement  of  the  tether  with  respect  to  the 
earth's  magnetic  field;  this  acts  to  drive  a  current  through  the  tether  to 
an  on-board  load.  Conversely,  if  current  is  forced  through  the  tether 
from  the  satellite,  the  electrodynamic  interaction  will  induce  a  force  and 
allow  the  tether  to  perform  as  a  so-called  Alfven  engine. 

The  concept  of  space  tethers  has  been  in  existence  for  about  a 
century,  with  the  earliest  reports  attributed  to  the  Russian  Konstantin  E. 
Tsiolkovsky  in  1895.  He  discussed  the  construction  of  a  tower  to  be  built 
at  the  equator  which  would  reach  to  geostationary  heights,  the  purpose  of 
which  was  to  create  an  environment  without  terrestrial  gravity.  The  next 
proposals  came  after  a  hiatus  of  about  sixty-five  years.  The  "heavenly 
funicular"  of  Y.N.  Artsutanov  and  the  "skyhook"  of  John  D.  Isaacs  each 
envisioned  the  tether  as  a  means  to  lift  payloads  into  earth  orbit. 
Additional  investigations  continued  throughout  the  sixties  and  seventies 
which  have  considered  the  concept  of  a  "space  necklace,"  essentially  an 
extension  of  the  heavenly  funicular  in  which  equally-spaced  funiculars 
about  the  equator  serve  as  payload-carrying  elevators  from  earth  to  GEO; 
and  a  "wheeled  tether",  in  which  long  rotating  tethers  are  placed  into  an 
orbit  such  that  at  specific  times  one  end  of  the  tether  would  just  touch 
the  earth's  surface  where,  since  it  would  briefly  appear  stationary 
relative  to  the  rotating  earth,  payloads  could  be  attached  to  be  slung  up 
into  orbit  as  the  tether  rotated  away. 

Initial  in-space  experiments  were  conducted  on  Gemini  XI  and  XII  in 
1966,  which  investigated  the  complex  dynamics  of  tethers  and  of  satellites 
attached  by  them.  The  first  practical  application  of  an  orbiting  tether 
in  LEO  was  developed  by  Mario  D.  Grossi  in  1972,  who  showed  that  the 
tether  could  be  used  as  an  antenna  to  generate  ultra  low  frequency  (ULF) 
emissions.  This  concept  was  further  developed  by  Giuseppe  Colombo  et  al 
in  1974  in  a  paper  which  proved  the  dynamic  feasibility  of  a  shuttle-borne 
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tether  system,  in  addition  to  several  other  additional  applications. 
Since  that  time,  several  other  endeavors  have  explored  the  problems  of 
tether  dynamics,  deployment  and  retrieval  of  long  tethers  from  a  free- 
flying  vehicle,  and  the  practicality  of  tethered  constellations  in  space. 

One  of  the  first  attempts  to  evaluate  the  electrodynamic  interaction 
of  long  conducting  tethers  and  the  ionosphere  was  a  paper  by  Dobrowolny, 
Colombo  and  Grossi  in  1976.  It  developed  computer  models  to  analyze  the 
distribution  of  the  induced  potential  along  the  tether  and  the  resulting 
current  from  this  potential.  Subsequently,  NASA  has  conducted  periodic 
seminars  and  sponsored  a  long  series  of  projects  to  investigate  the  many 
and  varied  facets  of  this  problem,  which  have  led  to  the  manifesting  of 
the  Tethered  Satellite  System  (TSS)  on  a  shuttle  flight  in  1991.  This 
experimental  package,  sponsored  by  the  U.S.  and  Italian  space  agencies,  is 
designed  to  evaluate  the  feasibility  of  generating  power  with  an 
electrodynamic  tether. 

The  purpose  of  this  thesis  is  three-fold:  first,  to  evaluate  the 
comparative  operation  of  an  electrodynamic  tether  which  is  perfectly 
insulated  from  the  surrounding  space  plasma  with  one  which  is  operated 
completely  bare;  second,  to  evaluate  the  effects  which  the  design 
parameters  of  the  tether  itself  and  of  the  orbit  in  which  it  is  placed 
have  on  the  efficiency  of  the  tether  system;  and  third,  to  evaluate  the 
effects  which  the  operation  of  the  tether  in  various  modes  have  on  the 
shape  of  its  orbit.  The  first  will  be  accomplished  by  performing  an 
analytical  evaluation  of  the  magnitude  of  current  lost  along  the  length  of 
an  uninsulated  tether  as  a  result  of  captured  positive  ions  and  secondary 
electron  emissions.  The  second  will  be  accomplished  as  a  numerical 
evaluation  of  the  thermal  balance  of  the  tether  system  by  computing  the 
temperature  of  the  tether  as  it  moves  through  its  orbit,  then  relating 
this  to  the  system  efficiency  by  means  of  the  fluctuations  in  conductor 
resistance,  which  is  a  function  of  temperature.  The  third  will  be 
accomplished  by  calculating  the  external  forces  which  act  on  the  tether, 
particularly  the  induced  electromotive  force  caused  by  the  tether's 
operation,  then  numerically  integrate  the  variational  equations  for  the 
orbital  elements  with  these  forces  present. 


2.  THEORETICAL  BACKGROUND 


2.1  INDUCED  VOLTAGE 

The  Lorentz  Relation  specifies  that  the  force  on  a  charge  q,  moving  at 
velocity  u  in  an  electric  field  E  and  magnetic  field  B  is  given  by 

F  =  q  (E  +  u  X  B)  [N/m2] 

To  evaluate  this  properly,  since  F,  B  and  q  are  all  independent  of  motion 
for  non-relativistic  velocities,  the  field  E  must  be  viewed  from  different 
frames  of  reference.  If  E',  then,  is  defined  as  the  field  when  viewed  in 
a  rest  frame,  and  E  is  defined  as  the  field  when  viewed  in  a  frame  which 
is  moving  at  y  relative  to  the  rest  frame,  this  indicates 

E'  +  uxB  =  E+  (u-y)  xB 


Thus, 


E  =  E'  +  y  X  B 

If  y  is  then  defined  as  the  mean  mass  velocity  of  the  plasma  surrounding  a 
rigid  wire,  and  since  in  the  plasma  frame  E'  must  be  small,  the  plasma 
will  experience  a  potential  across  the  length  of  the  tether  of  (y  x  B)  * 
L.  This  may  be  expressed  as 

Voc  =  (v  X  B)  •  L  [V] 

where  L  represents  the  length  of  the  wire  and  y  is  the  velocity  of  the 
tether  relative  to  the  plasma.  If  no  current  flow  is  permitted  due  to  the 
presence  of  a  large  impedance  somewhere  in  the  circuit,  a  potential  of 
this  magnitude  will  develop  across  that  impedance  due  to  the  velocity  and 
magnetic  field  interaction.  This  impedance  could  be  at  the  points  of 
plasma-wire  contact,  in  the  wire  itself,  or  in  a  localized  electrical  load 
or  open  switch;  any  of  these  impedances,  or  a  combination  of  them,  will 
support  the  y  ^  B  •  L  voltage. 


Additionally,  the  force  induced  by  a  current  carrier  moved  through  a 
purely  magnetic  field  is 

F  =  (I  X  B)  L  [N] 

where  L  is  the  tether  length,  and  the  direction  of  F  is  dependent  on  the 
direction  of  I.  Thus,  if  current  is  allowed  to  pass  through  a  conducting 
tether  which  is  being  moved  through  a  magnetic  field  (e.g.  the  earth's) 
such  that  a  force  is  induced  which  is  opposite  to  the  direction  of  motion, 
this  force  is  termed  electromagnetic  drag;  if  the  induced  force  is  in  the 
direction  of  motion,  it  is  termed  electromagnetic  thrust.  Its  effect  is 
to  "steal"  or  "add"  energy,  respectively,  from  the  satellite's  orbit  and 
alter  its  orbital  altitude. 

Power  may  be  generated,  or  thrust  produced,  by  capitalizing  on  these 
aspects  of  the  Lorentz  Relation,  provided  that  current  can  be  caused  to 
flow  through  the  tether.  A  current  loop  may  be  established  by  "clamping" 
the  two  ends  of  the  tether  to  the  plasma  so  that  the  voltage  drop  between 
each  respective  end  and  the  plasma  potential  is  minimized.  As  the  plasma 
by  its  nature  is  highly  conductive,  the  "clamping"  will  effectively  create 
a  current  loop.  Theoretical  work  indicates  that  this  clamping  is  possible 
(although  with  varying  degrees  of  resistivity)  and  that  the  ionosphere 
will  pass  this  current  with  only  negligible  resistance.  Once  this  loop 
has  been  established,  the  mechanical  power  added  to  the  system  may  be 
expressed  as 

P  =  F  •  V  =  L  (I  X  B)  •  V 

=  I  (L  X  B)  •  V  =  -I  (y  X  B)  •  L 

which  can  be  seen  to  be  the  product  of  the  current  and  the  open  circuit 
voltage.  This  expression,  integrated  over  time,  represents  the  energy 
which  can  be  extracted  from,  or  injected  into,  the  orbit  of  the  tether. 
Of  course,  only  a  fraction  of  this  can  be  delivered  to  a  load;  conversely, 
a  driving  power  source  would  have  to  provide  more  than  this  amount. 


2.2  EQUIVALENT  CLOSURE  CIRCUIT 

In  order  to  produce  a  current  from  the  potential  difference  created  by 
the  electromotive  force,  the  ends  of  the  tether  must  make  electrical 
contact  with  the  earth's  plasma  environment.  Once  this  contact  is 
established,  the  current  flows  through  a  'phantom'  loop  consisting  of  the 
tether  itself  and  the  ionospheric  plasma.  The  medium  in  which  the  tether 
is  orbiting  may  be  modelled  as  a  set  of  parallel  transmission  lines,  the 
magnetic  flux  tubes,  which  are  touched  sequentially  by  the  ends  of  the 
tether  as  it  moves  in  its  orbit.  At  this  touch,  electrons  leave  the 
positive  end  of  the  conducting  tether  and  enter  the  ionospheric  plasma 
where  they  produce  charge  packets  which  travel  along  these  magnetic  flux 
tubes  as  a  plasma  wave  in  a  variety  of  modes,  mostly  as  low-frequency 
Alfven  waves.  The  charge  packets  from  the  opposite  tether  ends  have,  of 
course,  opposite  polarities,  and  they  must  eventually  neutralize  each 
other.  However,  this  is  of  no  concern  to  the  spacecraft,  since  the 
distance  travelled  by  the  plasma  wave  down  each  magnetic  tube  in  the  time 
the  tether  contactor  is  in  touch  with  it  is  on  the  order  of  the  contactor 
size  times  the  ratio  of  the  wave  speed  to  the  contactor  speed,  and  is  no 
more  than  a  few  hundred  meters  in  any  case.  The  nature  of  the  eventual 
"closure"  in  the  E-layer  of  the  ionosphere  is  mainly  of  interest  in  the 
sense  of  providing  a  diagnostic  tool  for  ground  observations  of  the 
phenomenon,  and  perhaps  as  a  concern  for  its  electromagnetic  noise 
emission.  The  ends  of  the  tether  act  as  sliding  contacts,  so  that  this 
phenomenon  in  the  ionosphere  is  fundamentally  AC.  However,  since  the 
ionosphere  is  a  continuous  medium,  the  ends  of  the  tethers  do  not  leave 
one  line  before  touching  the  next  one,  so  that  the  current  in  the  tether 
itself  is  DC  and  nearly  constant. 

The  basic  equation  of  the  circuit  is 


Vqc  =  IR  +  Vioad  '/anode  ^cathode  '/ion 


where 


Vqc  =  Induced  voltage  across  the  tether  =  v  x  B  •  L 
I  =  tether  current 
R  =  resistance  of  the  tether 
Vload  =  voltage  drop  across  the  load 

Vanode  "  voltage  drop  across  the  anodic  contactor 

Vcathode  =  voltage  drop  across  the  cathodic  contactor 

Vi  on  =  voltage  drop  through  the  ionosphere 

Additionally,  there  are  current  leakage  losses  along  the  length  of  the 
tether  which  are  dependent  on  the  degree  of  effectiveness  of  the 
insulation  of  the  conducting  wire,  and  the  inductance  which  would  arise 
from  the  part  of  the  current  carrying  tether  which  remains  wound  around  a 
deployer  reel.  A  schematic  of  this  circuit  path  is  shown  in  Figure  1. 
Obviously,  if  either  the  conducting  wire,  the  ionosphere  or  the  contact  of 
the  tether  and  plasma  do  not  readily  pass  the  current,  i.e.  if  the 
effective  impedance  of  any  of  these  is  high,  then  the  efficiency  of  the 
system  will  plummet.  It  thus  becomes  essential  to  minimize  those  voltage 
drops  in  the  system  besides  the  load. 

In  performing  the  analysis  of  this  system,  the  impedances  of  this 
circuit,  with  the  exception  of  the  impedance  due  to  the  reeled  tether 
inductance,  were  modelled  mathematically  to  produce  a  simplified 
equivalent  circuit.  It  is  assumed  that  the  tether  is  completely  deployed, 
so  that  the  impedance  due  to  inductance  is  zero.  A  discussion  follows  of 
the  various  models  used  and  assumptions  made  to  obtain  this  equivalent 
ci rcuit. 


2.2.1  Ionospheric  Impedance 


The  power  losses  in  the  ionosphere  can  be  properly  called  'radiation 
losses',  and  have  been  evaluated  in  detail  in  References  15  and  18.  These 
analyses  made  calculations  of  radiation  resistance  in  the  framework  of 
linear  cold  plasma  theory,  referring  in  particular  to  low  frequency  Alfven 
waves  and  to  waves  in  the  vicinity  of  the  ion  cyclotron  frequency.  The 
conclusion  of  Reference  18  was  that  the  effective  radiation  power  loss  for 
a  tether  of  length  20  km  would  be  minimal,  and  implied  that  a  maximum 
total  radiation  resistance  of  approximately  of  1.2  ohms  was  to  be 


expected.  Reference  18  also  provided,  after  a  much  more  rigorous  analysis 
based  on  antenna  theory,  an  alternate  means  of  determining  the  approximate 
radiation  impedance  in  the  Alfven  regime.  The  author  here  showed  that  the 

impedance  of  an  orbiting  long  wire  is  about  where  n^  is  the  Alfven 
refractive  index.  This  index  varies  from  day  to  night,  typically  falling 
in  the  range  of  100  by  night  to  1000  by  day,  indicating  an  impedance  of 
about  1  ohm  by  night  and  .1  ohm  by  day,  and  thus  corroborating  the  first 
estimation. 

Other  analyses  have  been  made  to  estimate  the  radiation  impedance 
caused  by  resonance  at  higher  frequencies  in  the  Whistler  range,  notably 
References  4  and  22.  Their  conclusions  are  that  impedances  on  the  order 
of  1000  ohms  should  be  expected  in  these  regimes,  and  that  the  critical 
factor  which  determines  the  magnitude  of  the  impedance  is  the  smallest 
dimension  of  the  wire.  Reference  22  concludes  that  this  is  the  diameter. 
Obviously,  impedances  of  this  magnitude  would  make  the  effectiveness  of  a 
tether  system  marginal  at  best.  More  recent  work,  as  detailed  in 
Reference  15,  makes  it  clear  that  the  critical  dimension  is  that  of  the 
current  interface,  i.e.,  the  contactor.  The  authors  of  References  4  and 
22  now  concur  on  this  point.  As  will  be  discussed  in  the  next  section, 
the  effective  size  of  the  plasma  contactors  is  generally  accepted  as  on 
the  order  of  50  meters,  which  would  indicate  an  impedance  of  about  1  ohm. 

The  conclusion  which  may  be  drawn  is  that  the  wave  impedance  of  the 
ionosphere  for  the  electrodynamic  tether  will  not  present  a  major  limiting 
factor  to  the  current  flow  through  the  circuit.  As  the  full  implication 
of  these  results  is  not  yet  well  understood  and  still  the  topic  of  some 
debate,  this  analysis  will  assume  this  conclusion  to  be  accurate  and 
presume  that  the  total  radiation  impedance  in  the  ionosphere  is 
consistently  less  than  10  ohms. 

2.2.2  Contactors 

Several  methods  of  'bridging'  the  plasma  sheath  at  the  ends  of  the 
tether  have  been  devised.  The  degree  of  success  in  effecting  good  contact 
between  the  tether  and  the  plasma,  i.e.  their  level  of  impedance,  vary 


markedly  among  the  several  methods.  What  follows  is  essentially  a  brief 
summary  of  these  different  plasma  bridges. 

Those  contactors  placed  at  the  positive  end  of  the  tether,  the  anode, 
are  essentially  electron  collectors  (or  ion  emitters)  and  may  be  either 
passive  or  active.  Conversely,  the  contactors  at  the  negative  end,  the 
cathode,  are  all  active  electron  emitters. 

Cathodic  contactors,  as  discussed  in  reference  10,  pose  little 
fundamental  difficulty,  since  electrons  can  always  be  injected  into  the 
ionosphere  with  little  energy  loss,  provided  the  emitter  doses  not  charge 
up  positively  in  the  process.  As  long  as  the  anodic  contactor  is  doing 
its  job  of  collecting  electrons  at  an  appropriate  rate,  this  charging  will 
not  occur.  The  hollow  cathode  is  the  device  almost  universally 
recommended  to  serve  as  the  plasma  bridge  at  the  cathodic  end  of  the 
tether.  This  device  injects  a  plasma,  typically  formed  from  a  neutral  gas 
such  as  mercury,  xenon,  or  argon,  into  the  ionosphere.  This  increases  the 
density  of  ions  in  the  immediate  vicinity  of  the  contactor,  thus  providing 
a  sufficient  thermal  electron  density  to  carry  the  full  tether  current  in 
either  direction  at  any  distance  from  the  tether  and  until  it  merges  into 
the  ambient  ionospheric  plasma  currents.  Being  active  plasma  contactors, 
hollow  cathodes  require  on-board  electrical  power  and  a  gas  supply  to 
operate.  Various  estimates  have  been  made  for  this  voltage  drop  and  the 
required  mass  consumption  rate,  notably  in  References  21  and  29.  These 
indicate  that  sufficient  current  levels  can  be  attained  with  voltage  drops 
on  the  order  of  100  volts.  The  attendant  mass  consumption  rate  can  be 
estimated  as 


where 

Ce  =  random  thermal  speed  of  electrons 
Vs  =  ambi polar  speed  of  sound 
a  =  degree  of  ionization 
m,  =  ion  mass 

This  yields  a  mass  flow  rate  to  current  ratio  on  the  order  of  2  kg/A/year 
for  10%  ionized  argon. 


The  primary  competitor  of  the  hollow  cathode  device  is  the  electron 
gun,  which  has  the  advantages  of  being  an  already  existing  technology 
which  utilizes  no  mass  flow  and  provides  a  ready  control  of  current,  but 
which  has  the  disqualifying  disadvantage  of  requiring  large  levels  of  on¬ 
board  electrical  power  to  draw  sufficient  levels  of  current,  due  to  space 
charge  limitations  at  the  gun  exit.  These  voltage  levels  have  been 
estimated  at  from  1  to  10  kilovolts  to  draw  20  amps  of  current,  which  is 
on  the  order  of  the  voltage  potential  induced  across  a  moderately  sized 
tether.  In  the  hollow  cathode,  this  space  charge  is  neutralized  by  the 
plasma  ions;  thus,  very  low  voltage  drops  are  possible. 

Anodic  contactors  may  either  collect  electrons  or  emit  ions.  The 
former  may  also  be  either  passive  or  active  collectors.  The  passive  type 
function  essentially  by  placing  a  surface  in  the  ionosphere,  the  size  of 
which  is  large  enough  to  collect  the  required  magnitude  of  current  from 
the  ambient  plasma.  Obviously,  if  the  ambient  plasma  density  is  low,  the 
size  of  the  surface  area  must  be  proportionally  increased.  As  an  example, 
if  the  ambient  plasma  density  is  loH  m"^  then  the  random  current  density 

(given  by  is  .8  mA/m^.  If  the  required  current  level  is  10 
amps,  this,  then,  dictates  a  surface  area  of  12500  m^.  Equivalently,  a 
surface  area  of  this  magnitude  represents  a  sphere  of  approximately  45 
meters  radius,  or  a  flat  screen  80  meters  square.  Diagrams  of  these 
configurations  are  shown  in  Figure  2.  It  is  apparent  that  if  a  thin- 
skinned  balloon  is  used  as  the  sphere,  it  will  be  both  heavy  and  draggy, 
since  it  will  present  a  large  surface  area  perpendicular  to  the  direction 
of  velocity.  If  a  mesh  of  thin  wires  is  used  as  the  flat  screen,  its  mass 
and  drag  can  both  be  substantially  decreased;  it  will  still,  however, 
present  a  dynamical  monstrosity.  As  an  alternative,  the  plasma  cloud 
released  by  a  hollow  cathode  device  can  be  used  as  a  substitute  for  this 
actual  surface  area.  This  would  clearly  improve  things  but  does  have  its 
own  limitations.  Reference  29  presents  calculations  which  indicate  that  a 
cloud  radius  of  about  40  meters  is  necessary  to  collect  10  amps  of 
current.  This  is  further  complicated  by  the  fact  that  electrons  must  be 
able  to  diffuse  across  the  magnetic  field  into  the  central  plasma 
contactor,  so  the  only  useful  portion  of  the  plasma  cloud  is  that  in  which 


the  collision  frequency  exceeds  the  gyro  frequency,  which  occurs  for 
either  high  plasma  densities  (~  3.5  x  10^7  n,-3)  and/or  high  plasma 
turbulence. 

The  primary  deficiency  in  all  of  these  types  of  electron  collectors  is 
that  they  are  totally  dependent  on  the  ambient  electron  density  in  order 
to  continually  provide  a  sufficient  current  level.  The  Smithsonian 
Astrophysical  Observatory  has  conducted  a  study  of  electron  density  versus 
orbital  time,  shown  as  Figure  3,  calculated  using  a  model  based  on  the 
three-dimensional  maps  constructed  by  Jones  and  Stewart  (Reference  25). 
This  model  is  limited  in  that  it  is  based  on  data  from  November  1966;  as 
this  date  is  in  the  middle  of  the  solar  cycle,  the  model  has  no  included 
seasonal  or  solar  activity  variations.  It  also  does  not  correctly  model 
the  density  variations  at  lower  altitudes,  but  is  accurate  enough  at  300 
km  to  show  the  problem  with  electron  collectors  as  anodic  plasma  bridges. 
Specifically,  the  ambient  electron  density  varies  by  as  much  as  two  orders 
of  magnitude,  from  10^2  ni~3  to  lO^^  m"3,  between  the  sunlight  and  shadow 
portions  of  the  system's  orbit.  At  the  low  points,  the  surface  area 
required  to  collect  sufficient  current,  either  physical  surface  or  that  of 
a  plasma  cloud,  becomes  huge. 

What  seems  to  be  a  solution  to  this  problem  is  to  replace  these 
electron  collectors  at  the  anode  with  an  ion  emitter.  The  obvious 
advantage  of  this  device  is  its  near  independence  from  the  ambient 
electron  density,  in  contrast  with  the  discussion  above.  Such  a  device 
would  be  essentially  an  ion  engine  minus  the  accelerator  grid,  which  would 
not  be  necessary  as  the  ions  would  leave  the  device  because  of  the 
potential  gradient.  The  mass  consumption  of  this  type  device  is  somewhat 
greater,  though  not  significantly,  than  that  of  the  hollow  cathode;  it  has 
been  estimated  in  Reference  29  as  approximately  2  kg/A/year  with  ammonia 
as  the  working  gas.  The  power  requirements  would  likewise  be  somewhat 
higher  than  the  cathode,  but  still  be  on  the  order  of  100  volts  to 
maintain  adequate  current  levels. 

For  this  analysis,  the  assumption  will  be  made  that  the  cathodic 
plasma  bridge  is  a  hollow  cathode  device,  and  that  the  anodic  bridge  is  a 
pure  ion  emitter.  Since  mathematical  models  for  the  impedances  of  these 


contactors  are  complex,  nonlinear,  and  unknown  functions  of  tether 
current,  it  will  be  assumed  that  the  voltage  drop  across  both  the  anode 
and  cathode  is  a  constant  100  volts  apiece.  Finally,  it  will  be  assumed 
that  each  device  has  sufficient  mass  in  place  for  unlimited  continuous 
operation  so  that  the  consideration  of  fuel  consumption  by  the  contactors 
may  be  neglected. 

2.2.3  Insulation 

In  order  to  operate,  the  conducting  material  of  the  tether  must  pass  a 
current,  the  magnitude  of  which  is  determined  by  the  amount  of  power  or 
thrust  to  be  produced.  If  the  conductor  is  covered  with  a  dielectric 
material  of  sufficient  thickness,  then  no  current  is  lost  to  the 
surrounding  plasma  along  the  length  of  the  tether.  If,  alternately,  the 
tether  is  left  bare,  then  the  current  is  reduced  by  an  amount  which  is  a 
function  of  the  tether's  length  and  diameter,  the  potential  bias,  and  the 
ambient  ion  density.  The  serious  disadvantage  of  an  insulated  tether  is 
that  most  dielectrics  tend  to  differentially  charge  due  to  ion  impacts, 
and  then  arc.  This  arcing  typically  occurs  at  potential  levels  well  below 
the  4  kilovolt  level  which  is  to  be  anticipated  for  a  20  kilometer  tether. 
Additionally,  if  pinholes  are  introduced  into  the  dielectric,  as  is 
expected  due  to  micrometeroid  impacts,  the  magnitude  of  the  current 
leakage  through  these  pinholes  is  on  the  same  order  of  the  magnitude  of 
current  lost  by  a  bare  tether.  Thus  it  appears  that  unless  a  dielectric 
can  be  chosen  which  is  resistant  to  this  arcing  and  which  is  durable 
enough  to  withstand  most  debris  impacts,  for  all  intents  and  purposes  the 
insulated  tether  will  act  electrically  like  a  bare  tether,  but  have  a 
greater  mass. 


Reference  18  indicates  that  the  necessary  thickness  of  polyethylene  to 
prevent  surface  discharge  is  .38  mm.  This  was  evidently  determined  in  a 
quick  and  dirty  way  by  dividing  the  maximum  potential  which  the  tether 
would  encounter  by  the  dielectric  strength  of  the  insulating  material  to 
obtain  the  required  thickness,  then  multiplying  by  an  arbitrary  safety 
factor.  The  thickness  indicated  may  be  obtained  using  values  of  4750 
volts  for  the  maximum  potential,  5  x  10^  V/cm  for  the  dielectric  strength 
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of  polyethylene,  and  a  safety  factor  of  4.  Reference  34  presents  a 
similar  analysis  to  determine  the  required  dielectric  thickness  to 
insulate  a  satellite  from  the  Jovian  flux. 

Obviously,  insulation  of  this  thickness  will  prevent  current  from 
leaking  directly  through  the  dielectric.  The  cost  of  this  is  primarily  in 
additional  mass  and  in  the  manufacturing  complexities  involved  in 
producing  such  a  tether.  For  an  aluminum  tether,  2  mm  in  diameter  and  20 
km  in  length,  using  polyethylene  as  the  insulator  will  add  approximately 
52.3  kg  to  each  tether  mass;  the  insulator  to  conductor  mass  ratio  is  thus 
.297.  Comparatively,  anodized  aluminum,  with  a  dielectric  strength  of  5  x 
10^  V/cm,  would  require  a  thickness  of  only  38  microns  to  insulate  the 
conductor,  and  consequently  an  insulator  to  conductor  mass  ratio  of  .041. 
These  mass  ratios  may  be  improved  somewhat  by  tapering  the  insulation  so 
that  it  is  thickest  at  those  points  on  the  tether  where  the  largest 
potential  gradient  is  expected  to  occur.  However,  as  stated  above, 
although  the  current  leakage  along  the  length  of  the  tether  would  be 
prevented  by  this  insulation,  the  major  problems  still  to  be  faced  would 
be  arcing  and  pinhole  leakage. 

Arcing  is  defined  in  Reference  17  as  the  rapid  rearrangement  of  charge 
by  punch-through  (internal  breakdown  of  dielectrics),  by  flash-over 
(dielectric  breakdown)  between  surfaces,  or  between  surfaces  and  space. 
The  mechanisms  leading  to  arcing  are  not  well  understood.  They  fall  into 
three  general  categories:  ungrounded  insulators;  buried  or  trapped  charge 
layers;  and  differential  voltage  distributions.  Typically,  the  charge 
build-up  associated  with  ungrounded  insulators  can  be  prevented  by  proper 
design  of  the  insulating  jacket.  Buried  charges  are  caused  by  mid  energy 
electrons,  on  the  order  of  50  KeV,  which  penetrate  the  insulating  jacket 
and  deposit  their  charge  in  the  dielectric.  Electron  energies  of  this 
order  are  not  typically  encountered  in  LEO,  the  operating  regime  of  the 
tether.  Differential  voltage  distributions  occur  either  when  the  surface 
of  the  dielectric  is  at  a  lower  potential  than  the  conductor  beneath  it, 
and  small  gaps  or  cracks  exist  in  the  dielectric  which  expose  the 
conductor  and  allow  charge  to  escape  locally,  or  through  charge  deposition 
on  the  surface  of  the  dielectric  due  to  the  ambient  plasma  and  due  to 


secondary  emissions  from  surface  impacts.  Protecting  the  tether  from  this 
phenomenon  requires  production  techniques  which  encase  the  tether  without 
any  gaps,  cracks,  or  crevices  in  the  dielectric,  in  addition  to  the 
selection  of  a  dielectric  which  has  a  high  dielectric  strength  and  is 
durable  enough  to  resist  erosion  and  penetration. 

The  consequences  of  a  bare  tether  will  be  discussed  in  a  subsequent 
section,  as  will  the  use  of  anodized  aluminum  as  a  dielectric  which 
appears  to  have  the  desired  characteristics  of  high  dielectric  strength 
and  durability,  in  addition  to  values  of  absorptance  and  emittance  which 
will  cause  the  tethe'  to  operate  in  a  favorable  temperature  range.  For 
the  remainder  to  this  analysis,  it  will  be  assumed  that  the  tether  is 
perfectly  insulated  by  a  jacket  of  anodized  aluminum,  and  that  arcing  and 
pinhole  damage  do  not  cause  any  difficulties  in  system  operation. 

2.2.4  Reduced  Circuit 

With  these  simplifying  models,  the  equivalent  circuit  to  be  analyzed 
is  as  shown  in  Figure  4.  Here,  the  constant  voltage  drops  across  the 
anode  and  cathode  have  been  subtracted  from  the  variable  emf,  Vqc*  The 
low  impedance  of  the  ionosphere  has  been  assumed  negligible,  and  therefore 
omitted.  The  tether  is  assumed  to  be  perfectly  insulated  so  that  there 
are  no  leakage  losses,  and  is  assumed  to  be  completely  deployed  so  that 
there  exists  no  reeled  inductance.  The  electromotive  force  remains 
variable  because  of  the  magnetic  field  variations.  The  tether  resistance 
also  is  variable  because  of  temperature  variations  throughout  the  orbit. 

2.3  MODEL  FOR  INDUCED  VOLTAGE 

Essential  to  the  evaluation  of  various  tether  configurations  is  the 
accurate  calculation  of  the  fluctuations  in  the  electromotive  force  (emf), 
defined  as  (y  x  B)  •  L.  If  it  is  assumed  that  the  orbit  of  the  tether  is 
circular,  v  becomes  a  constant  fixed  by  the  orbital  altitude.  Further,  if 
variations  in  L  are  computed  as  a  consequence  of  temperature  changes,  the 
single  variable  to  be  evaluated  becomes  B. 


The  magnetic  field  data  used  in  the  analyses  of  the  thermal  variations 
which  the  tether  experiences,  which  were  all  conducted  for  an  invariant 
circular  orbit,  were  produced  by  the  SAO  using  the  NEWMAG  software  code. 
This  Is  a  spherical  harmonic  model  of  the  earth's  magnetic  field  which 
starts  with  the  assumption  that  the  predominant  portion  of  the  earth's 
magnetic  field  B  can  be  expressed  as  the  gradient  of  a  scalar  potential 
function  V  as 


B  =  -VV 


and  that  the  expression  of  the  potential  V  of  the  internal  field  can  be 
written  as  a  series  of  spherical  harmonics: 


^  nil  Jo  K  *  "n 


where  a  is  the  equatorial  radius  of  the  earth;  g  and  h  are  Gaussian 
coefficients;  r,  6,  and  ^  are  the  geocentric  distance,  coelevation,  and 
east  longitude  from  Greenwich,  respectively,  which  define  a  point  in 
space;  and  are  the  associated  Legendre  functions.  The  Gaussian 
coefficients  are  determined  empirically  by  a  least  squares  fit  of 
measurements  of  the  magnetic  field.  Together  with  the  associated  Legendre 
functions,  it  is  thus  possible  to  calculate  the  magnetic  field  at  any 
point  in  space. 


From  V,  three  orthogonal  components  of  the  magnetic  field  may  be  derived 
as : 


„  -1  av  ;  fa1nt2  P  f.ti.m _ .  .  ^n.m  _,1  ap"'"'(e) 

-  as  '  -nEl  '7'  raEo  19  cos  *  h  •  sin  m^J 

n.  -1  0V  -1  "  faln+2  P  _r  _n,m  _ ^  .  L.n.m  „n.m,«x 

“I*  •  TTiiS  a^  ■  iTiS  n5l  ra  m5o  "'1-9  *  P  '  cos  m^J  P  (6) 


The  data  produced  by  SAO  with  this  code  utilized  the  set  of  coefficients 
of  the  International  Geomagnetic  Reference  Field  (IGRF)  for  epoch  1975.0. 
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In  addition  to  the  determination  of  the  Gaussian  coefficients  as  mentioned 
above,  the  coefficients  of  the  I6RF  assume  that  the  associated  Legendre 
functions  are  Schmidt  normalized,  i.e.: 


[p”(e)]  sine  ds  .  -L.^. 

where  (5}j^  is  the  Kronecker  delta.  These  coefficients  together  with  the 
geocentric  system  of  reference  (coelevation  ff,  east  longitude  and 
geocentric  distance  r)  then  provide  the  vectorial  components  Bg,  B^,  and 
at  the  specified  point  where  the  electrodynamic  tether  is  located  in 
space. 


It  should  be  noted  that  this  model  approximates  the  various  magnetic 
anomalies  observable  in  LEO  which  would  not  be  accounted  for  had  a  pure 
tilted  dipole  model  been  used.  The  tilted  dipole  is  essentially  the 
spherical  harmonic  field  model  expanded  to  the  first  degree  (n  =  1)  and 
all  orders  (m  =  0,1).  As  can  be  seen  from  the  plots  of  (y  x  B)  •  L  versus 
time  for  a  24  hour  period  in  Figures  5A  and  6,  the  output  using  the  tilted 
dipole  model  is  simply  a  sinusoidal  oscillation  specified  by  the  cross 
product  of  y  and  B  (assuming  the  tether  is  perfectly  rigid  along  the 
radial  line),  where 

g^'^  sin5  -  [g^'^  cos^  +  h^'^  sin^j  cos0  j 

Here  is  the  equatorial  radius  of  the  earth,  r  is  the  orbital  radius,  g 
and  h  are  the  IGRF  coefficients  for  the  first  degree  and  orders  0  and  1,  0 
=  sin"l  (sin  i  sin  nt),  and  ^  =  sin"l  (cot  i  tan  nt),  where  i  is  the 
tether's  orbital  inclination  to  the  magnetic  polar  axis,  which  varies 
during  the  course  of  a  day,  and  nt  is  the  true  anomaly  of  the  tether  in 
its  orbit.  This  sinusoidal  variation  is  due  to  the  rotation  of  the  dipole 
axis  about  the  earth's  geographic  axis  during  the  course  of  a  day;  if  the 
earth's  rotation  were  neglected,  the  induced  voltage  potential  would 
remain  constant  throughout  the  orbit.  By  comparison,  the  output  from  the 
NEWMAG  code  shows  this  same  daily  sinusoid  due  to  the  earth's  rotation, 
but  additionally  very  large  intra-orbital  fluctuations  superimposed  on 


this  sinusoid  as  a  result  of  flying  over  or  near  the  actual  magnetic 
anomalies  in  the  earth's  field,  notably  the  great  negative  South  Atlantic 
anomaly  and  the  positive  Australian  and  Asian  anomalies.  For  the  analyses 
conducted  to  determine  the  variation  in  the  orbital  elements,  an  expansion 
to  the  third  order,  the  so-called  octopole  model,  was  used.  It  can  be 
seen  from  Figure  5B  that  this  represents  a  reasonably  close  approximation 
of  the  fully  expanded  model,  while  requiring  less  computational  time  and 
permitting  deviations  from  a  circular  orbit. 

The  inputs  for  NEWMAG  were  generated  by  the  Smithsonian  Astrophysical 
Observatory  (SAO)  using  their  SKYHOOK  computer  program.  This  software 
package  is  designed  to  model  the  electrodynamic  interactions  of  a  tethered 
satellite  system,  and  includes  models  for  the  earth's  magnetic  field  as 
described  above,  electron  density,  and  charge  collection  by  the  wire  and 
various  contactors.  It  functions  in  three  steps.  The  first  is  the 
establishment  of  initial  conditions  for  the  computer  run.  The  second  is 
the  analytical  calculation  of  orbital  motion  to  obtain  the  behavior  of  the 
tether  as  a  function  of  time.  The  third  is  the  analysis  of  the  output 
using  various  post-processors.  For  this  analysis,  a  rigid  tether  was 
assumed  and  its  mass  modeled  by  two  end  points.  The  initial  conditions 
for  the  position  and  velocity  of  each  of  these  mass  points  were  obtained 
using  a  pre-processor  called  DUMBEL,  which  takes  orbital  parameters  and 
computes  a  set  of  initial  conditions  which  allow  the  formulation  to  be 
started  in  tension  equilibrium  so  that  there  are  no  initial  oscillations 
in  the  radial,  in-plane,  or  out-of-plane  directions.  The  computations  by 
SKYHOOK  produce  the  state  vector  of  the  system  at  each  time  step,  which 
includes  the  position  and  velocity  of  each  mass  point  in  geocentric 
coordinates;  this  state  vector  is  the  input  for  NEWMAG. 

A  post-processor  written  by  SAO  then  took  the  SKYHOOK  values  for  the 
state  vector  produced  by  the  integration  to  calculate  the  vector  product 
Vqc,  ^nd  the  NEWMAG  values  for  the  magnetic  field  components,  as  described 
above,  to  calculate  the  magnitude  of  the  magnetic  field  strength,  at  50 
second  time  intervals.  These  two  data  were  the  basis  for  the  thermal 
analysis  conducted  for  varying  tether  surface  conditions  and  operating 
modes. 


3.  METHODS  OF  ANALYSIS 


In  this  examination  of  thermal  and  certain  dynamical  effects  on  the 
operation  of  an  electrodynamic  tether  system,  three  areas  in  particular 
will  be  di scussed. 

First,  an  analysis  of  the  consequences  of  leaving  the  conducting 
material  of  the  tether  bare,  instead  of  insulated  as  is  generally  presumed 
in  most  analyses,  will  be  made.  This  will  be  done  to  a  first 
approximation  by  calculating  the  magnitude  of  the  current  which  could  be 
expected  to  be  lost  along  the  length  of  an  uninsulated  tether  due  to  the 
capture  of  ions  from  the  surrounding  plasma  and  additionally  to  the 
emission  of  secondary  electrons  as  a  result  of  these  captures.  The  ratio 
of  this  current  loss  to  the  system  current  without  losses  will  provide  an 
effective  efficiency  for  an  uninsulated  tether. 

Second,  using  the  reduced  closure  circuit  described  in  Section  2.4 
above  as  the  basis  for  the  analysis,  an  examination  of  the  effects  of  the 
variation  of  the  operational  temperature  on  tether  performance,  as  defined 
by  its  efficiency,  as  a  result  of  the  revolution  of  the  tether  about  the 
earth,  will  be  made  by  numerically  integrating  the  temperature  equation 
forward  in  time  from  an  arbitrarily  selected  set  of  initial  conditions. 
The  orbital  elements  of  the  tether’s  orbit  will  be  varied  to  determine 
their  individual  effects,  as  will  the  design  parameters  of  the  tether 
itself.  This  analysis  will  be  restricted  by  the  assumptions  of  a  rigid 
tether  moving  in  a  perfectly  circular,  invariant  orbit.  Additionally,  the 
tether  will  be  simulated  to  be  operating  as  a  power  generator,  a  thruster, 
and  various  combinations  of  the  two. 

Third,  the  dynamical  effects  of  the  various  modes  of  tether  operation 
on  the  orbit  of  the  tether  will  be  examined  by  computing  those  external 
forces  which  act  on  the  tether  and  serve  as  perturbations  to  the  orbital 
elements.  Again,  these  perturbation  equations  will  be  numerically 
integrated  forward  in  time  to  determine  if  any  particular  method  of 
operation  is  physically  prohibitive. 
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3.1  INSULATION 


As  discussed  in  Section  2.2.3  above,  the  conducting  material  of  the 
tether  may  be  insulated  from  the  surrounding  plasma  to  eliminate  any 
current  leaking  along  its  length  by  coating  it  with  a  sufficient  thickness 
of  some  type  of  dielectric.  However,  a  tether  insulated  in  this  manner 
experiences  the  disadvantages  of  arcing  due  to  differential  charging  of 
the  dielectric,  and  large  leakage  currents  through  pinholes  in  the 
dielectric.  An  alternative  to  this  problem  is  to  actually  leave  the 
tether  bare  if  the  amount  of  current  which  leaks  along  its  length  is  of  a 
small  enough  magnitude  to  constitute  an  acceptable  loss.  The  following  is 
an  analysis  to  estimate  the  magnitude  of  this  current  leakage  loss. 

An  upwardly  deployed  tether  moving  in  a  prograde  orbit  will  charge 
negatively  with  respect  to  the  surrounding  plasma.  Thus,  the  current  loss 
which  the  tether  will  experience  along  its  length  will  be  the  result  of 
positive  ions  from  the  plasma  which  are  captured  by  the  negatively  biased 
tether,  plus  the  secondary  electrons  emitted  by  the  wire  upon  ion  impact. 
The  number  of  ions  which  are  captured  is  a  function  of  the  ion  density  at 
the  operational  altitude,  the  bias  of  the  tether,  the  ion  mass  and 
velocity,  and  the  thickness  of  the  tether. 

There  are  several  mathematical  models  which  may  be  used  to  estimate 
the  magnitude  of  the  current  collected  by  the  tether  as  it  moves  through 
the  surrounding  plasma.  The  most  appropriate  for  our  case  is  the  so- 
called  inertia  limited  current  collection  model,  in  which  current 
collection  is  due  only  to  those  ions  whose  trajectories  under  the 
attraction  of  the  wire  intercept  the  wire  itself.  It  acts  as  an  upper 
limit  for  the  amount  of  current  which  can  be  collected,  because  the 
presence  of  earth's  magnetic  field  restricts  somewhat  the  number  of  ions 
which  can  actually  strike  the  tether.  It  is  also  less  than  that  of  a 
space-charge  limited  collection  model,  in  which  a  positive  ion  sheath  is 
formed  around  the  tether  such  that  the  only  current  collected  is  due  to 
those  ions  which  randomly  "walk"  across  the  sheath  boundary.  As  detailed 
in  Reference  17,  this  type  of  mode  is  generally  more  applicable  to 
structures  which  have  larger  relative  dimensions  than  a  tether. 


If  sheath  effects  about  the  tether  are  ignored,  which  may  exaggerate 
the  estimate  of  ions  captured  as  described  above,  this  number  of  ions  may 
be  calculated  by  balancing  the  momentum  of  the  ion  far  from  the  tether 
with  the  momentum  of  one  which  just  grazes  the  tether's  surface.  As  shown 
in  Figure  7,  if  x  is  the  distance  from  a  line  passing  through  the  tether's 
center  which  is  parallel  to  the  velocity  direction,  and  the  velocity  of 
the  ion  far  from  the  tether  is  represented  as  Vq,  then 

xvq  =  rvi 

where  r  is  the  radius  of  the  tether  and  v^  is  the  velocity  of  the  ion  as 
it  grazes  the  tether.  The  momentum  balance  then  yields 

miVo^  =  mi  Vi^  -  2  e  V 

where  V  represents  the  magnitude  of  the  tether  bias.  Substituting  the 
first  equation  into  the  second  and  solving  for  x  shows  the  distance  from 
the  tether  greater  than  which  the  ions  will  not  impact,  and  hence  not 
cause  a  current  loss,  when  they  are  attracted  by  the  negative  bias.  Thus 


r  2 
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Obviously,  this  analysis  assumes  that  the  plasma  potential  far  from  the 
tether  is  zero.  This  is  a  reasonably  good  assumption  so  that  the  estimate 
for  X  should  also  be  good.  Having  determined  this  distance,  the  current 
loss  due  to  ion  capture  is  estimated  as 

IlOSS  =  2  X  L  ni  Vo  e 

where  L  is  the  length  of  the  tether  and  ni  is  the  ion  density  at  the 
operational  altitude.  Having  assumed  a  completely  neutral  plasma,  i.e.  n^- 
=  ng,  then 
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I,  =  2  e  r  L  n^  V 


Finally,  since 


>  >  1,  then 


2eV 


where  d  is  the  tether  diameter.  Using  typical  values  of  ne  =  10^1  m"3,  d 
=  3  mm,  L  =  20  km,  V  =  4  kV,  and  =  .016  kg/mol,  this  would  indicate  a 

current  leakage  of  approximately  .21  amps  if  the  tether  potential  with 
respect  to  the  surrounding  plasma  were  constant  along  its  length.  For  a 
tether  operating  to  generate  10  kw  of  power,  this  represents  an 
approximately  7.5%  loss  above  any  other  system  losses.  However,  since  the 
tether  potential  with  respect  to  the  plasma  really  varies  along  its 
length,  a  better  estimate  for  the  loss  can  be  obtained  by  defining  this 
variation  for  a  specified  system. 

For  the  case  of  a  tether  which  is  being  operated  to  generate  power, 
the  current  in  the  system  is  positive  upward,  as  shown  in  Figure  9,  where 
9A  shows  the  generator  mode  and  9B  the  thruster  mode.  From  9A  it  can  be 
seen  that  the  voltage  induced  across  the  length  of  the  tether,  Vqci  is 
equal  to  the  sum  of  the  potential  drops  across  the  contactors  at  the  anode 
and  cathode,  the  load  and  the  tether  itself.  This  assumes,  of  course, 
that  these  are  the  only  system  losses  as  discussed  in  the  formulation  of 
the  equivalent  circuit  above.  If  it  is  assumed  that  the  potential  of  the 
tether,  is  linear  with  z,  then 

=  Az  +  B 

Solving  for  the  integration  constants,  we  note  from  Figure  9A  that  the 
boundary  conditions  may  be  written  as 

at  z  =  0:  =  ^Po  -  -  Vl  =  B 

and 

at  z  =  L:  =  ^Po  -  Vqc-  +  Vc  =  Az  +  B 
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which  yields 


A  =  (Vl  +  Va  +  Vc  -  Voc)  /  L 

Therefore 

-  (V^  *  (f)  +  -  ''a  -  \ 

Using  this,  and  defining  the  plasma  potential  as  function  of  z  as 

h  =  ho  -  VcB  z 

where  ^po  is  the  plasma  potential  at  z  =0,  Vq  is  the  orbital  velocity, 
and  B  is  the  magnetic  flux,  we  can  then  write  an  equation  for  the 

potential  of  the  tether  with  respect  to  the  plasma  at  any  point  as 

=  ''l  *  "a  -  (''l  "  "a  *  ''c)  (!)  <2> 

Rewriting  Equation  1  for  an  arbitrary  value  of  z,  we  get  a  current  per 

unit  length  as 

^  ^  "e  d  j 

If  the  voltage  drop  due  to  the  tether's  resistance  is  written  as 


01 


V  '^cJ 


Here,  rj  = 


a  ''c  Vi 

''a  ''c 

I  -  V  I  3nd  S  =  ^ 


S[U6^-5^}-[U6^l2 


Examining  the  two  terms  of  this  equation,  it  can  be  seen  that  the 
first  is  equivalent  to  that  of  the  current  flow  due  to  the  voltage  drop 

across  the  tether  itself,  I  =  ^  ,  while  the  second  is  equivalent  to  that 
of  the  current  flow  due  to  the  ion  capture.  The  difference  in  sign 
between  the  two  terms  indicates  that  the  ion  capture  is  effectively  a 
current  loss. 

The  ratio  of  these  two  terms  is  essentially  an  efficiency  for  the  bare 
tether.  A  good  approximation  for  this  is  made  by  expanding  the  second 
term  as  a  binomial  series  and  neglecting  all  factors  of  order  greater  than 
one  and  the  products  of  all  factors.  This  neglect  is  reasonable  as  6a  ~ 
6r  «  1.  This  expansion  yields  a  ratio  expressed  as 
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Calculating  this  loss  using  the  typical  values  stated  previously,  and 
assuming  a  system  efficiency  of  90%,  a  value  of  3.24%  is  attained.  This 
represents  the  effective  current  loss  due  solely  to  ion  capture  along  the 
length  of  the  bare  tether,  and  is  about  half  of  what  was  calculated 
previously  for  a  tether  of  uniform  potential  along  its  length. 


The  loss  calculated  above  represents  only  that  current  which  is  due  to 
ion  collection.  The  total  current  loss  is  the  sum  of  several  additional 
components,  specifically  secondary  electron  emission  due  to  ion  impacts, 
photoelectric  emission,  backscattering,  and  thermal -induced  currents.  In 
a  low  earth  orbit,  such  as  that  in  which  the  tether  would  operate,  these 
components  are  typically  negligible  with  the  exception  of  secondary 
electron  emission.  So,  to  complete  this  estimate  for  current  leakage  from 
an  uninsulated  tether,  the  additional  current  losses  resulting  from 
secondary  electron  emission  due  to  ion  impacts  must  be  accounted  for. 


Since  the  emission  of  an  electron  is  the  equivalent  of  the  capture  of 
a  position  ion,  these  secondary  electrons  also  represent  an  effective 
leakage  loss.  References  11  and  17  both  discuss  this  and  present  graphs 
of  electron  emission  from  nickel  and  aluminum  targets.  For  ion  impacts  at 
the  4  kV  energy  level,  the  ratio  of  the  electron  emission  current  to  the 
positive  ion  current  is  approximately  60%  from  the  nickel  target,  and  240% 
from  the  aluminum  target.  Additionally,  as  corroboration,  it  is  indicated 
that  the  number  of  secondary  electrons  emitted  per  primary  impact  from  an 
aluminum  target  is  about  4  times  the  number  emitted  from  a  nickel  target. 
Hence,  the  total  current  leakage  loss  would  be  3.4  times  the  loss 
indicated  previously.  Together  with  the  ion  current  calculated  above,  a 
total  current  leakage  loss  due  to  positive  ion  capture  and  secondary 
electron  emissions  amounts  to  approximately  11.0%. 

Finally,  before  deciding  if  a  loss  of  this  magnitude  is  acceptable 
with  respect  to  the  mass  savings,  an  estimate  should  be  made  to  determine 
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the  degree  of  conductor  loss  over  time  due  to  the  ion  impacts.  The  rate 
of  mass  loss  may  be  expressed  as 

^  -  k'  *^1  ^i 

dt  I  e 

where  is  the  leakage  current  as  determined  above,  and  K  is  a  constant, 
characteristic  to  the  conductor,  representing  the  number  of  conductor 
atoms  ejected  per  ion  impact.  Data  from  Reference  11  may  be  extrapolated 
to  make  an  estimate  for  K  for  an  aluminum  target  as  approximately  .44  A1 
atoms  /  O'*"  ions.  This  is  probably  a  pessimistic  estimate  as  an  oxide 
coating  has  been  found  to  develop  on  bare  aluminum  in  the  presence  of 
atomic  oxygen,  such  as  in  the  ionosphere,  which  sputters  very  little. 
Thus,  the  degree  of  mass  loss  over  the  expected  operational  life  of  the 
tether  is  about  2  kg  in  10  years.  Hence,  conductor  mass  loss,  or 

sputtering,  does  not  appear  to  pose  a  problem  for  an  uninsulated  tether. 


A  consistent  analysis  may  be  made  for  a  tether  operated  in  the 
thruster  mode  by  referring  to  the  potential  diagram  in  Figure  9B.  The 
result  of  this  process  is  an  estimate  for  the  leakage  loss  due  to  positive 
ion  capture  as 


Loss 


where  now  V  represents  the  voltage  drop  across  the  source,  5=^,5= 
w  y  s  a  V  c 

c  oc 

,  ri  =  —  .  The  percent  loss  when  operated  in  this  mode,  including 
losses  for^secondary  electron  emissions,  amounts  to  approximately  13.6%. 


Similarly,  an  analysis  can  be  made  for  a  downwardly  deployed,  and 
hence  positively  biased,  tether.  The  expressions  for  the  leakage  losses 
are  similar  to  those  derived  above.  Now,  however,  the  capture  of  the 
negative  electrons  in  the  surrounding  plasma  are  the  cause  of  the  leakage, 
and  so  the  mass  of  the  electron  must  be  substituted  into  the  expressions 
in  place  of  the  ion  mass.  Since  the  ratio  of  the  mass  of  atomic  oxygen  to 
that  of  an  electron  is  approximately  30000,  this  drives  the  loss  ratio  to 
100%.  A  positively  biased,  uninsulated  tether  is  thus  prohibitive. 


As  a  consequence  of  this  analysis,  an  effective  efficiency  for  the 
uninsulated  tether  may  be  defined  as 

’?Loss  =  1  -  3.4  Loss 

where  Loss  is  the  ratio  calculated  above  to  estimate  the  current  loss  due 
to  ion  capture  along  the  length  of  the  tether  compared  to  the  current  flow 
without  losses,  and  the  3.4  accounts  for  the  additional  losses  due  to 
secondary  electron  emission.  As  previously  stated,  however,  for  the 
remainder  of  this  analysis  it  will  be  assumed  that  the  tether  is  perfectly 
insulated  so  that  no  current  leakage  losses  occur. 

3.2  METHOD  OF  THERMAL  ANALYSIS 

Given  the  reduced  circuit  as  discussed  in  Section  2.4  above  so  that 
the  equation  of  the  circuit  is  written  as 

IR  =  Voc  -  200  -  ^  [V] 

and  conducting  the  system  analysis  for  a  constant  load  power,  Pl,  the  only 
remaining  variables  to  be  determined  are  the  tether  resistance,  R,  and  the 
system  current,  I.  The  resistance  of  the  tether  is  primarily  a  function 
of  the  temperature  of  the  conductor;  hence,  a  thermal  analysis  will  be 
done  to  calculate  the  temperature,  and  subsequently  the  resistance, 
electromotive  force,  and  efficiency,  as  a  function  of  orbital  time.  The 
current  will  be  determined  such  that  a  constant  load  power  is  maintained 
at  an  arbitrarily  specified  level  when  the  tether  is  operating  as  a  power 
generator.  When  the  tether  is  operated  as  a  thruster,  the  current  will  be 
determined  by  the  thrusting  mode,  i.e.,  to  maintain  a  constant  current 
level  while  thrusting,  to  maintain  a  constant  force  level  while  thrusting, 
or  to  maintain  a  constant  voltage  level  while  thrusting. 

In  the  analysis,  the  parameters  of  the  tether  itself,  i.e.,  the 
length,  diameter,  type  of  conductor  material,  and  surface  characteristics, 
will  be  examined  to  determine  their  relative  effects  on  the  system 
operation.  Likewise,  the  orbital  parameters,  i.e.,  the  altitude, 
inclination,  distance  from  the  sun,  and  orientation  to  the  sun,  will  be 
examined  to  determine  their  relative  effects. 
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3.2.1  Thermal  Balance 


As  in  any  system,  the  thermal  balance  of  the  tether  may  be  expressed 


^in  ^  "^gen  ’  ^out  ~ 

Formally,  gout  consists  of  radiation  to  free  space,  radiation  to  other 
objects,  and  conduction  losses;  pgen  is  the  ohmic  heat  generated  as 
current  passes  through  the  tether;  and  q^-p  consists  of  solar  irradiation, 
earth  irradiation,  earth  albedo,  self-irradiation,  and  aerodynamic 
heating.  For  this  analysis  it  will  be  assumed  that  Qout  consists  only  of 
radiation  to  free  space,  and  that  self-irradiation  is  a  negligible  portion 
of  qip.  These  assumptions  are  fairly  well  justified.  The  neglect  of 
self-irradiation  follows  from  the  assumption,  made  previously,  of  a 
perfectly  rigid  tether.  Radiation  from  the  tether  to  free  space  dominates 
that  to  other  nearby  objects,  primarily  the  base  station,  again  as  a 
consequence  of  the  rigid  tether  assumption;  thus,  only  considering 
radiation  to  free  space  seems  reasonable.  Neglecting  conduction  losses  is 
not  as  easily  justified.  However,  the  ends  of  the  tethers  must  be 
electrically  insulated  to  prevent  current  from  discharging  through  the 
base  supports.  As  electrical  insulators  are  generally  also  good  thermal 
insulators,  conduction  between  the  tether  and  supports  should  not  occur. 

Given  these  assumptions,  the  components  of  the  thermal  balance  are  as 
fol low: 


A  r  ^ 

q  -  m  C  dt 


ear,  AT  =  ear,  ir  d  L  T 
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"solar  ■  S  "  L  Z 


^aero  "  2  ^  '^rel 

'’albedo  =  -28  S  T  d  L  F  Icos  ?l 

'’earth  =  237  e  ir  d  L  F 
These  equations  use  variables  defined  as: 

d  =  tether  dian-eter 
L  =  tether  length 
T  =  temperature 
=  Boltzman's  constant 
I  =  current 

Os  =  solar  absorptivity 
e  -  emissivity 
Z  =  sun  factor 
F  =  earth  view  factor 
Vf-el  =  relative  velocity 
/7a  =  atmospheric  density  at  altitude 

?  =  biased  true  anamoly 
p  =  resistivity 

t  =  time 

m  =  conductor  mass 

C  =  specific  heat 

S  =  solar  constant  =:  1368  W/m2 

More  detailed  explanations  of  the  selections  of  values  for  these 
parameters  will  be  discussed  in  Section  3.2.2. 

The  thermal  capacitance  of  the  system  must  be  examined  to  determine 
how  quickly  the  tether  will  reach  thermal  equilibrium,  and  thus  to  decide 
whether  the  equilibrated  state  can  be  used  as  a  point  for  evaluation. 
Writing  the  equation  above  for  the  thermal  balance  strictly  in  terms  of 
temperature,  we  have 


where  A’  =  e  ir  d  L,  and  B'  =  q-jn  +  Qgen*  Assuming  that  B'  is  a 
constant  with  respect  to  temperature,  this  implies  that  the  time  constant 
for  the  system  is 


T  =  - - 3 

4  e  Jr  d  L 

where  Tq  is  the  equilibrium  temperature.  Using  typical  values  for  these 
parameters,  as  follow: 


m  =  400  kg 
C  =  875  J/kg/OK 
d  =  3  mm 
L  =  20  km 
To  =  240  OK 
£  =  .85 

this  yields  a  time  constant  of  t  =  700  seconds,  or  about  11.6  minutes. 
Compared  to  an  approximately  90  minute  orbit,  this  indicates  that  a  state 
of  equilibrium  would  be  inappropriate  for  evaluation,  and  that  the  thermal 
transients  caused  by  the  capacitance  of  the  tether  must  be  taken  into 
account. 

Despite  this  transience  throughout  the  orbit,  at  any  given  time  the 
tether  may  be  assumed  to  be  at  a  uniform  temperature.  Since  in  our  case, 
the  length  of  the  tether  will  always  be  much  longer  than  the  radius,  the 
tether  may  be  considered  a  cylinder  with  a  uniformly  distributed  heat 
source  (the  current).  Thus,  it  will  be  assumed  throughout  this  analysis 
that  all  portions  of  the  tether  are  simultaneously  at  the  same  temperature 
when  performing  computations. 

3.2.2  Calculation  of  Components  of  Thermal  Balance 

Having  made  these  assumptions,  an  examination  of  the  separate 
components  of  the  thermal  balance  can  be  conducted. 


Radiation  from  the  tether  to  free  space  is  solely  a  function  of  the 
temperature  of  the  conductor,  as  specified  by  the  Stefan-Bol tzmann  law, 
presuming  that  the  other  variables  in  the  equation  may  be  expressed  in 


terms  of  T.  While  this  is  indeed  the  case,  the  variations  of  the 
emissivity  and  the  diameter  are  so  slight  over  the  range  of  temperature 
variations  which  the  tether  will  undergo  that  they  may  be  considered  as 
constants.  The  length,  on  the  other  hand,  may  vary  significantly  over  an 
orbit  depending  on  the  magnitude  of  the  temperature  variation  and  must  be 
taken  into  account.  This  variation  in  length  may  be  expressed  as 

T  =  «e  (^2  -  h) 

where  Oq  is  the  coefficient  of  expansion,  typically  about  12.9  x  10'^ 
m/m/OK  for  aluminum. 

Solar  irradiation  of  the  tether  is  a  function  of  the  intensity  of  the 
solar  flux  and  the  magnitude  of  the  tether's  surface  area  presented 
normally  to  the  incident  radiation.  The  solar  flux  is  a  parameter  which 
is  typically  considered  to  be  a  constant;  it  does,  however,  vary  somewhat 
within  the  11-year  solar  cycle  and  routinely  is  affected  by  solar  storms. 
The  average  value  of  the  solar  flux  outside  of  the  earth's  atmosphere  over 
the  past  several  years  has  been  about  1368  W/m^.  Since  the  solar  flux 
varies  as  the  inverse  square  to  the  distance  from  the  sun,  there  is  a 
variation  of  ±3?^  due  to  the  eccentricity  of  the  earth's  orbit  throughout 
the  year.  The  maximum  occurs  at  the  perihelion  of  the  earth's  orbit 
(about  3  January,  1399  W/m^);  the  minimum  at  aphelion  (about  4  July,  1309 
W/m^) .  The  absorptance  of  the  material  is,  like  the  emittance,  relatively 
invariant  with  such  temperature  variations  as  the  tether  will  experience, 
and  will  also  be  a  considered  a  constant.  What  here  has  been  called  the 
sun  factor  is  obtained  by  a  series  of  rotation  angles  to  calculate  the 
tether  area  perpendicular  to  the  incoming  solar  rays  at  a  given  point  in 
the  orbit.  The  angles,  as  shown  in  Figure  10,  represent  the  inclination 
of  the  orbit,  i,  relative  to  the  equatorial  plane;  the  obliquity  of  the 
equatorial  plane,  p,  relative  to  the  ecliptic  plane;  the  longitude  of  the 
line  of  notes,  [I,  relative  to  the  normal  to  the  earth-sun  line;  and  the 
true  anomaly  of  the  tether  in  its  orbit,  f,  measured  from  the  line  of 
nodes . 


Both  the  inputs  due  to  earth  irradiation  and  albedo  involve  the  earth 
view  factor  for  the  tether.  This  view  factor  is  a  function  of  the 
altitude  of  the  tether  in  its  orbit,  and  its  orientation  with  respect  to 
the  earth-tether  radial  line.  Formulae  for  various  orientations  are  given 
in  Reference  10.  Assuming  that  the  tether  remains  perfectly  aligned  with 
this  radial  line,  the  view  factor  may  be  calculated  as 

P  ^  Q  -  sin  6  cos  0 

t 

where  6  =  sin 
At  an  altitude  of  300  km,  this  yields  an  earth  view  factor  of  F=. 31413. 

The  earth  albedo  flux  is  essentially  that  fraction  of  the  solar 
irradiation  reflected  back  onto  the  tether  as  a  result  of  scattering  in 
the  atmosphere  and  reflection  from  clouds  and  earth  surfaces.  This 
fraction  is,  much  like  the  solar  flux  itself,  basically  a  constant, 
although  variations  do  occur  as  a  result  of  changes  in  cloud  cover  and 
snow  and  ice  covered  regions.  Typical  values  in  the  low  latitudes  are  30% 
±2%.  The  primary  difference  between  the  two  is  that  instead  of  only  the 
normal  area  being  irradiated,  here  the  entire  surface  area  of  the  tether 
receives  the  incident  albedo,  mitigated  by  the  earth  view  factor  and  the 
proximity  of  the  tether  to  the  terminator.  This  proximity  is  included  as 
the  cosine  of  the  biased  anomaly,  such  that  the  albedo  is  a  maximum  when 
directly  facing  the  sun  and  decreasing  to  zero  at  the  point  in  the  orbit 
where  the  earth's  shadow  is  entered. 

Irradiation  by  the  earth  itself  involves  essentially  an  "earth  flux", 
caused  by  solar  radiation  which  has  been  absorbed  by  the  earth  and 
atmosphere  being  reemitted  as  thermal  radiation.  Again,  this  remains 
approximately  constant  at  237  ±7  W/m^,  with  the  earth  view  factor 
mitigating  the  total  surface  area.  Now,  however,  the  emissivity  of  the 
tether  surface  determines  the  energy  input  instead  of  the  absorptance  as 
in  the  solar  irradiation  and  albedo  calculations.  This  is  because  the 
earth's  radiation  is  not  at  the  wavelengths  of  the  solar  radiation;  hence, 
the  solar  absorptance  cannot  be  used  in  this  calculation. 


Aerodynamic  heating  is  a  function  of  the  density  of  the  atmosphere  at 
the  operational  altitude  (which  in  turn  is  a  function  of  atmospheric 
temperature  and  the  solar  cycle);  and  the  velocity  of  the  atmosphere 
relative  to  the  tether.  Each  of  these  tend  to  vary  as  the  satellite  moves 
in  its  orbit.  In  performing  the  calculation  for  drag  the  maximum  values 
which  presented  themselves  were  used  in  order  to  produce  a  conservative 
estimate. 

Empirical  formulas  have  been  developed  for  the  calculation  of  density. 
For  an  orbital  altitude  of  300  km  this  is: 
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For  this  evaluation,  the  exosphere  temperature  was  selected  to  be  1600“  K, 
the  maximum  for  the  solar  cycle.  Thereafter,  in  subsequent  calculations 

=  6.225  x  10"^^  kg/m^ 

was  used. 

Atmospheric  velocity  varies  according  to  latitude  and  altitude  as 
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where  f  is  the  true  anomaly  and  the  28.5“  is  the  assumed  orbital 
inclination  used  in  this  analysis.  The  velocity  of  the  tether  relative  to 
the  surrounding  air  is  then  a  vector  difference,  which  may  be  expressed  as 


2  2  2 

V  i=v  .+  v.-2v  uV.  cosrf 
rel  orb  air  orb  air  " 


where  ^  may  be  expressed  as  28.5“  sin  f.  For  our  purposes  then,  the 
maximum  relative  velocity  at  an  orbital  altitude  of  300  km  is 


Vel  =  i 
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Using  these  values  for  and  v^eli  the  maximum  energy  input  from 
aerodynamic  heating  is  given  by 

qaero  =  12-108  d  L  [W] 

Obviously,  when  calculated  in  this  manner,  the  value  for  aerodynamic 
heating  will  vary  throughout  the  orbit. 

The  ohmic  heat  generated  by  the  current  passing  through  the  tether 
varies  with  both  temperature  and  with  the  magnitude  of  the  open  circuit 
voltage,  (v  x  B)  •  L  or  Vqc-  Since  the  electrical  resistance  may  be 
expressed  as 

R  =  ^ 

r  d 

where  both  p  and  L  are  functions  of  temperature,  the  resistivity  as 

/3  =  (l  -  0-^  (t  -  293.15°]] 

and  L  as  described  previously,  and  since  I  varies  throughout  the  orbit 
proportionally  with  Vqc.  this  Joule  heating  will  be  a  constantly  varying 
source  of  energy  for  the  tether.  The  magnitude  of  the  current  will 
further  depend  on  the  amount  of  power  which  is  required  to  be  delivered  to 
the  load. 

Using  these  equations  and  substituting  into  the  equation  for  the 
thermal  balance,  the  temperature  transients  that  would  occur  for  tethers 
of  various  diameters  as  the  system  orbits  the  earth  and  the  changes  in 
those  parameters,  such  as  resistance,  efficiency,  and  current,  which  are 
functions  of  temperature  can  be  evaluated.  In  this  analysis,  the 
efficiency  of  the  tether  itself,  not  including  other  system  losses,  may  be 
expressed  as 


As  can  be  seen  from  Figures  11,  which  represent  the  variations  in  the 
properties  of  a  tether  operating  as  a  pure  generator,  temperature 
variations  of  up  to  40  degrees  can  be  expected  throughout  the  orbit, 
depending  on  the  diameter  of  the  tether  under  consideration.  Similarly, 
the  efficiency  can  fluctuate  by  about  30  points.  Efficiency  improves  when 
the  tether  is  operated  at  lower  temperatures  (an  expected  result  as 
resistivity  increases  with  temperature),  and  when  the  current  levels  are 
low  so  that  the  ohmic  losses  are  minimized. 

3.2.3  Conductor  Material 

The  choice  of  the  conducting  material  of  the  tether  must  be  made  on 
the  basis  of  the  material's  conductivity  and  density.  A  lighter  tether  is 
preferable  to  a  heavier  one  from  the  standpoint  of  the  cost  of  energy  to 
place  the  tether  into  orbit;  however,  a  light  tether  which  is  highly 
resistive  will  have  to  be  made  thicker  in  order  to  pass  a  specified  level 
of  current.  Conversely,  a  more  conductive  tether  will  decrease  the  ohmic 
losses  in  the  system,  hence  improving  the  system’s  overall  efficiency; 
however,  if  it  is  an  extremely  dense  material  the  improvement  in 
efficiency  will  not  compensate  for  the  increased  energy  costs. 

To  this  end,  the  parameter  to  evaluate  is  the  ratio  of  conductivity 
and  density,  ^  ,  such  that  the  determination  of  merit  is  for  a  large 
value.  Examinations  of  this  parameter  were  made  for  various  materials,  as 
indicated  in  the  table  below,  which  were  selected  on  the  basis  of  either 
high  conductivity,  low  density,  or  both.  Both  density  and  conductivity 
were  allowed  to  vary  over  the  range  of  temperatures  in  which  the  tether 
would  most  probably  operate,  although  the  variations  in  density  of  any  of 
the  materials  examined  over  the  small  range  of  probable  temperatures  was 
negligible.  The  characceristics  of  the  material  evaluated  are  also  listed 
in  the  table.  The  variations  in  the  parameter  are  shown  graphically  for 
the  several  materials  in  Figure  12. 


Conductor  Characteristics 


Type 

Ref  Resistivity 

Coeff  Resistivity 

Ref  Density 

Coeff  Density 

Aluminum 

2.857 

x 

10-8 

.00424 

2707 

-6.912  X  10-5 

Copper 

1.724 

x 

10-8 

.0068 

8910 

0 

Si  1 ver 

1.629 

x 

10-8 

.0061 

7148 

0 

Nickel 

7.8 

x 

10-8 

.0069 

8890 

0 

Platinum 

10. 

X 

10-8 

.00393 

2145 

0 

It  can 

be  seen 

from 

Figure  12  that 

the  parameter  for  aluminum  is 

significantly  better  (higher)  than  those  of  the  other  materials  evaluated. 
Thus,  aluminum  will  be  the  conducting  material  used  in  this  analysis  of  a 
tether  system. 

3.3  ORBITAL  VARIATIONS 

When  an  external  set  of  accelerations,  in  addition  to  the  primary 
acceleration  due  to  the  main  gravitational  field,  acts  on  an  orbiting 
satellite  the  parameters  which  define  the  orbit  vary  with  time.  Appendix 
C  details  the  theoretical  development  which  leads  up  to  the  variational 
equations  of  Lagrange  and  Gauss  which  may  be  used  to  calculate  the  rate  of 
change  of  the  orbital  elements.  If  the  external  forces  acting  on  the 
satellite,  in  this  case  the  tether,  can  be  calculated,  then  the  vector 
components  of  that  total  force  will  specify  how  the  orbit  will  be  altered 
in  the  presence  of  that  force. 

3.3.1  Disturbing  Acceleration  Vector 

The  disturbing  acceleration  is  essentially  the  sum  of  all  external 
forces  which  act  on  the  tether  system  throughout  its  orbit.  The 
components  which  were  included  for  this  analysis  were  the  induced 
electromagnetic  force,  aerodynamic  forces,  solar  radiation  forces,  tether 
1 ibration-induced  forces,  and  the  effects  of  the  earth's  oblateness. 
Additionally,  there  are  external  forces  caused  by  induced  eddy  currents  in 
the  system  which  may  interact  with  the  earth's  magnetic  field,  and  by  the 
solar  wind  and  cosmic  dust;  however,  these  are  negligible  by  comparison 
with  the  other  forces  listed  here  and  will  not  be  included. 


For  this  analysis,  an  earth-centered  coordinate  system  was  used  with 
the  first  coordinate  axis  parallel  with  the  eccentricity  vector,  the  third 
parallel  to  the  angular  momentum  vector,  or  normal  to  the  orbital  plane, 
and  the  second  such  that  it  formed  a  right-hand  set  with  the  other  two, 
and  consequently  is  parallel  to  the  parameter  of  the  orbit.  Although  this 
is  apparently  a  rotating  coordinate  frame  with  respect  to  inertial  space, 
the  calculation  of  the  variational  equations  as  described  in  Appendix  C  is 
made  instantaneously  along  an  osculating  orbit,  so  that  the  Coriolis, 
Euler,  and  centripetal  accelerations  are  not  present.  This  coordinate 
system  is  depicted  in  Figure  13.  Thus,  we  can  write 


cos  f 

V  cos  f  -  V.  sin  f 

r  =  r 

sin  f 

and  y  = 

v[^  sin  f  +  v^  cos  f 

0 

0  ‘ 

with  v^  representing  the  radial  component  of  velocity,  equal  to  ^ , 
and  v^  representing  the  component  normal  to  v^  in  the  orbital  plane,  equal 
fn  e  cos  f) 

The  disturbing  acceleration  due  to  the  induced  electromagnetic  force 
may  be  expressed  as 

ad  =  g  X  B)L 
m 

where  L  is  the  tether  length,  m  the  system  mass,  I  the  vector  represen¬ 
tation  of  the  current,  and  B  the  magnetic  field  strength.  For  this 
analysis,  it  was  assumed  that  the  system  could  be  modeled  as  a  point  mass; 
hence,  the  tether  is  rigid  and  oscillates  only  in  the  orbital  plane. 
Further,  it  was  assumed  that  a  positive  I  represents  current  flowing  up 
the  tether,  i.e.  the  generator  mode.  The  oscillation  of  the  tether  was 
modelled  in  accordance  with  Reference  26,  assuming  that  the  length  of  the 
tether  remains  constant  and  that  the  moment  of  inertia  of  the  tether  about 
the  longitudinal  axis  is  negligible  compared  to  the  moments  of  inertia 
about  the  two  transverse  axes.  Assuming  that  e  «  1,  as  would  be  the  case 
for  the  orbit  of  the  tether.  Equation  5.102b  of  this  reference  becomes 
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where  a  represents  the  in-plane  libration  angle;  the  particular  solution 
of  this  equation  is  a  =  e  sin  nt.  Hence,  approximating  f  =  nt,  we  can 
define  the  current  vector  as 

cos  a  cos  f  -  sin  a  sin  f 
I  =  I  cos  a  sin  f  +  sin  a  cos  f 

0 

in  terms  of  the  e,  p,  h  coordinate  system.  The  magnetic  field,  B,  was 
modeled  for  this  analysis  by  an  octopole  model,  in  which  the  field  model 
is  expanded  to  the  third  degree  and  all  orders,  with  the  vector  components 
again  referred  to  the  e,  p,  h  coordinate  system  described  above.  This  is 
essentially  as  described  in  Section  2.3;  a  representation  of  the  field 
strength  computed  by  this  model  appears  as  Figure  5B.  A  comparison  with 
Figure  5A,  the  simple  tilted  dipole  model,  and  Figure  6,  the  fully 
expanded  model,  shows  that  the  octopole  is  significantly  better 
approximation  than  the  dipole  which  does  not  require  the  computational 
effort  of  the  expanded  model. 

Typically,  the  magnitude  of  the  electromagnetic  acceleration,  for  a 
current  level  of  approximately  1  amp  and  a  tether  length  of  20  km,  is  on 
the  order  of  1  N/kg. 

The  component  of  the  disturbing  acceleration  due  to  aerodynamic  forces 
was  calculated  as  a  function  of  the  density  of  the  atmosphere  at  the 
operational  altitude  (which  in  turn  is  a  function  of  atmospheric 
temperature  and  the  solar  cycle);  the  velocity  of  the  atmosphere  relative 
to  the  tether;  and  the  drag  coefficient  for  the  system.  Values  for 
aerodynamic  heating  were  calculated  in  the  formulation  of  the  thermal 
balance;  values  for  aerodynamic  drag  may  be  made  using  the  same 
parameters,  as  follows. 

Like  electromagnetic  drag,  the  aerodynamic  drag  is  a  force  which  is 
distributed  evenly  along  the  tether's  length.  It  may  be  calculated  by  the 
i ntegral 
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Assuming,  as  we  have,  that  all  of  the  variables  in  the  integral  remain 
constant  over  the  length  of  the  tether,  the  drag  force  may  be  expressed  as 

td  -  i  S  '’a  4e1  “  L 

where  d  is  again  the  tether  diameter  in  meters,  and  Fj  represents  the  drag 
on  the  tether  acting  opposite  the  direction  of  the  relative  velocity.  For 
a  tether  diameter  of  three  millimeters,  the  magnitude  of  this  force  is  on 
the  order  of  .02  Nt. 

Vectorially,  Reference  10  defines  the  relative  velocity  as 

vrel  =  Y  -  Yair 

where  v  is  the  tether's  orbital  velocity,  as  defined  above,  and  ya-|>  is 
the  local  air  velocity.  This  local  air  velocity  is,  in  turn,  defined  as 

=  [fI  ^ 

ej 

where  X  is  the  local  latitude,  and  Vq  is  the  rotational  velocity  of  the 
earth  at  the  equator.  These  calculations  are  made  in  the  program  attached 
as  Appendix  B,  referring  the  vector  components  again  to  the  e,  p,  h 
system.  The  air  density  is  calculated  as  a  function  of  both  altitude,  H, 
and  exosphere  temperature,  Tex#  formulae  listed  in  Reference  10. 

These  are 


for  70  <  H  <  118 


[H  -  95) 
2600 

[l.47  X  10 


0  T^^  pOOO  -  T^J 

2.9  (H  -  200)110 

T  .. 


for  118  <  H  <  200 


for  H  >  200 
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with  Tex  varying  in  relation  to  the  tether's  position  with  respect  to  the 
sun  as  Tgxo  (1  *12  sin  f ) .  With  the  exception  of  H,  X,  Texo»  Vorbi 

r,  all  other  variables  were  assumed  constant  for  this  calculation.  Those 
not  specified  by  the  dimensions  of  the  tether  were  chosen  as  Cq  =  2.2;  and 
Tex  =  1100  K. 

The  component  of  the  disturbing  acceleration  due  to  the  solar 
radiation  force  is  given  by 

Fr  =  -K  A  P  iy 

where  K  is  a  dimensionless  constant  which  indicates  the  reflectivity  of 
the  material;  A  is  the  satellite's  area  perpendicular  to  the  sun;  and  P  is 
the  solar  momentum  flux,  approximately  4.4  x  10“0  kg/m/s^.  As  indicated, 
the  force  acts  opposite  the  direction  of  the  sun.  Typically,  K  varies 
between  0  and  2,  with  a  value  of  0  indicating  a  translucent  material,  and 
2  indicating  a  totally  reflective  material.  Using  a  dark  coating,  such  as 
anodized  aluminum  as  will  be  discussed  in  Section  4.1.1,  would  set  K  at 
about  1.25.  As  before,  A  would  equal  d  L  and  must  be  referred  to  the  x, 
y,  z  coordinate  frame  to  be  consistent.  Thus,  the  greatest  magnitude  of 
this  force  becomes  =  .0005  N. 

To  evaluate  the  effect  of  the  earth's  oblateness  and  transform  it  into 
a  disturbing  acceleration,  it  is  necessary  to  express  it  first  as  a 
disturbing  function,  then  find  its  gradient.  If  the  distribution  of  mass 
in  the  earth  is  assumed  to  be  symmetric  about  the  polar  axis,  so  that  all 
tesseral  harmonics  are  neglected,  and  only  the  effects  of  the  first  three 
zonal  harmonics  of  the  earth's  gravitational  potential  are  included,  then 
the  disturbing  function  is 

I*  =  -(?]  nh  \  (r)"  *'> 

where 

cos  ^  •  1^ 

The  assumption  of  axial  symmetry  is  reasonably  good  for  the  earth  and  is 
routinely  used;  likewise,  the  approximation  of  the  earth's  oblateness  by 
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only  the  J2,  J3,  and  J4  terms  is  also  routinely  used  and  is  good  since  the 
coefficients  for  the  earth  become  progressively  smaller  (J3  is 
approximately  500  times  smaller  than  J2). 

The  calculation  of  the  gradient  of  R  with  respect  to  ri  and  hence  the 
disturbing  acceleration,  is  detailed  in  Appendix  B.  The  magnitude  of  this 
acceleration  is  typically  about  .01  N/kg. 

Due  to  the  gravity  gradient  which  occurs  in  extended  bodies,  such  as  a 
tether,  small  net  forces  are  produced  as  the  body  librates  during  its 
orbit.  As  detailed  in  Reference  26,  and  assuming  for  this  analysis  that 
the  tether  is  represented  as  a  simple  system  of  two  equal  masses,  these 
forces  may  be  written  as  a  radial  force 


-1.5/i  m  r,  ,  r  .  2  ) 

^  -  4 -  U  -  1.5  sin  a) 

r 

and  an  azimuthal  force 

where  a  is  the  libration  angle  in  the  plane  of  the  orbit,  equal  to  e  sin 
nt  as  described  above.  To  first  order,  then,  the  disturbing  acceleration 
due  to  tether  libration  may  be  written  as 


2 

cos  f  +  e  sin  f 
sin  f  (1  -  e  cos  f) 
0 


in  the  e,  p,  h  frame.  The  magnitude  of  this  acceleration  is  thus 
typically  on  the  order  of  10'^  N/kg. 


The  sum  of  these  individual  disturbances  is  the  disturbing 
acceleration  vector  to  be  used  in  the  calculation  of  the  variation  of  the 
orbital  parameters. 
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The  equations  which  were  used  to  analyze  the  effects  of  the  disturbing 
acceleration  on  the  orbit  of  the  tether  were  the  variational  equations  in 
vector  form,  as  are  described  in  Appendix  C.  These  are 


da  ,  2^ 
dt 


^  ^  r  sin  6  . 
dt  h  sin  i  -h 


di  r  cos  6 


ih  •  ^d 


ad  -  cos  1  ^ 


^1-  +  ^  a 
^2  3y  ^d 


de  1 


dt  =  ;  [2(!;  •  5d)  I  -  (I  •  5d)  “-(£•!!)  id] 


where  the  perturbative  derivative  of  the  true  anomaly  is 


)[( 


h  (cos  f  +  el 
P 


and  the  scalar  derivative  of  the  eccentricity  is 

^  _  (I  •  »]  [r  ■  id]  *  [p"  -  •  id) 

dt  /*  a  e 

Although  these  equations  are  valid  for  any  coordinate  frame,  as  long  as  r, 
V,  and  aj  are  consistent,  it  is  obvious  that  a  numerical  solution  using 
them  will  blow  up  when  applied  to  orbits  characterized  by  very  low 
inclination  and/or  eccentricity.  Of  course  this  is  because  the  line  of 
nodes,  and  hence  (1,  is  undefined  for  zero  inclination,  or  equatorial , 
orbits;  similarly,  u  is  undefined  for  zero  eccentricity,  or  circular, 
orbits. 


4.  APPLICATION  OF  METHODS  &  INTERPRETATION  OF  RESULTS 


Each  of  the  methods  of  analysis  described  in  Section  3  were  applied  to 
a  tether  system  using  the  basic  computer  programs  included  in  Appendices  A 
and  B  with  parameters  varied  to  achieve  desired  conditions  of  operation. 
The  system  was  modelled  to  maintain  a  constant  power  to  the  load  in  all 
modes  of  operation.  The  reactions  of  the  tether,  as  evidenced  by 
variations  in  operating  temperature  and  in  orbital  elements,  when  operated 
under  varying  conditions  were  then  examined  to  determine  if  any  of  these 
operating  conditions  posed  a  problem  which  would  make  the  operation  of  the 
tether,  with  that  condition  present,  prohibitive. 

4.1  THERMAL  EFFECTS 


Using  the  method  of  thermal  analysis  described  in  Section  3.2,  the 
relative  effects  of  certain  design  parameters  of  the  tether  system  may  be 
examined  to  determine  the  degree  of  their  influence  on  the  system's 
operation.  Generally,  these  design  parameters  fall  into  two  groups; 
first,  those  related  to  the  physical  tether  itself,  such  as  the  surface 
characteristics,  conductor  material,  diameter  and  length;  and  second, 
those  related  to  the  type  of  orbit  in  which  the  tether  is  operating,  such 
as  the  altitude,  inclination,  longitude  of  the  line  of  nodes,  and  relative 
position  with  regard  to  the  sun. 

This  examination  will  briefly  describe  the  relative  effects  of  these 
parameters  on  the  average  operating  temperature  of  the  system,  and 
graphically  present  the  variations  in  system  temperature,  current  level, 
tether  resistance,  magnitude  of  induced  electromagnetic  force,  and  system 
efficiency. 

Prior  to  this,  however,  it  is  useful  to  examine  how  the  properties  of 
a  typical  tether  vary  throughout  several  orbits.  As  described  in  Section 
3.2.1,  far  from  being  in  a  state  of  equilibrium  the  tether's  temperature 
is  continually  changing  as  it  orbits  the  earth.  In  a  typical  orbit,  the 
variation  in  its  operating  temperature  will  be  as  depicted  in  Figure  15A. 
As  can  be  seen,  the  peak  temperatures  occur  when  the  greatest  surface  area 


is  exposed  to  the  sun.  This  occurs  in  two  places,  just  before  and  just 
after  'noon'.  Also,  as  shown  in  Figure  15B,  tether  resistance  follows 
this  temperature  variation.  Thus,  it  is  expected  that  if  current  is  held 
constant,  the  system's  efficiency  will  also  follow  this  variational 
pattern.  This  can  be  seen  in  Figure  15C.  If  current  is  allowed  to 
fluctuate,  as  in  the  case  of  operation  as  a  pure  generator  which  is 
required  to  maintain  a  constant  load  power,  this  fluctuation  will 
basically  follow  the  fluctuations  in  the  induced  voltage,  shown  previously 
in  Figure  6.  In  this  case,  the  efficiency  tends  to  follow  the  current 
fluctuations. 

4.1.1  Effects  of  Tether  Parameters 
Absorptance  &  Emittance 

As  can  be  seen  from  the  equation  listed  in  Section  3.2.2,  the 
characteristics  of  the  surface  of  the  tether  will  primarily  determine  the 
operating  temperature  of  the  system. 

Of  the  external  flux  which  impinges  on  the  tether,  only  aerodynamic 
heating  is  not  a  function  of  either  absorptance  or  emittance.  Thus,  we 
would  expect  that  if  the  tether  surface  was  characterized  by  high 
emittance  and  high  absorptance,  then  a  large  amount  of  the  incoming  flux 
would  be  absorbed  and  the  tether  would  heat  up.  Conversely,  the  tether 
sheds  energy  by  radiating  it  to  free  space:  as  shown  before,  this 
radiation  is  proportional  to  the  surface  emittance.  Thus,  a  tether 
characterized  by  a  high  emittance  would  rid  itself  of  heat  more  easily 
than  one  with  low  emittance.  We  would  thus  expect  that  the  greatest 
operating  temperatures  would  appear  for  tethers  which  possess  a  high 
absorptance,  to  absorb  energy,  but  low  emittance,  to  ineffectively  emit 
energy.  Similarly,  the  lowest  average  operating  temperatures  should  occur 
for  tethers  characterized  by  high  emittance  and  low  absorptance. 

Figure  16A  shows  the  average  operating  temperatures  for  identical 
tethers  with  varying  surface  characteristics.  As  expected,  those  tethers 
which  operate  at  the  lowest  temperatures  are  characterized  by  values  of 
low  absorptance  and  high  emittance,  while  those  operating  at  the  highest 
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temperatures  are  those  characterized  by  high  absorptance  and  low 
emittance.  The  operating  temperature  levels  are  most  influenced  by  the 
absorptance  of  the  material. 

It  is  convenient  to  use  the  trends  of  Figure  16A  to  require  the  tether 
to  operate  in  a  specific  temperature  range  by  choosing  a  surface  material 
which  possesses  the  correct  characteristics  to  achieve  the  desired 
temperature  level.  As  stated  previously,  the  tether  will  operate  more  and 
more  efficiently  as  the  ohmic  losses  are  reduced.  Since  the  system  power 
is  a  function  of  the  current  level,  which  thus  cannot  be  reduced  without 
also  reducing  this  power,  this  is  accomplished  by  minimizing  the 
resistance  of  the  conductor.  The  lowest  resistances  will  occur  at  low 
temperatures,  as  resistivity  is  linearly  proportional  with  temperature. 
Thus,  it  is  desirable  to  cause  the  tether  to  operate  at  as  low  a 
temperature  as  possible.  Choosing  a  surface  material  characterized  by 
high  emittance  and  low  absorptance  will  accomplish  this.  Figure  16B 
graphically  shows  the  characteristics  of  some  typical  dielectrics  as 
determined  by  their  values  of  absorptance  and  emittance.  As  can  be  seen, 
a  good  choice  for  the  tether  surface  would  be  anodized  aluminum;  this  will 
be  the  assumed  surface  in  the  remainder  of  the  analysis. 

Conductor  Material 

As  discussed  in  Section  3.2.3,  the  selection  of  the  conductor  material 
for  the  tether  is  determined  primarily  by  a  conductance  to  mass 
comparison.  Figure  12  shows  this  parameter  for  some  typical  conducting 
materials.  Figure  17  compares  the  average  operating  temperature  of 
similarly  dimensioned  tether  made  of  copper.  As  stated  before,  aluminum 
represents  the  best  choice  for  the  anticipated  range  of  tether  operating 
temperatures,  in  addition  to  causing  the  tether  to  operate  at  a  low 
temperature. 

Tether  Diameter 

With  the  exception  of  the  ohmic  heat  generated  internally  by  the 
tether,  the  variation  in  the  average  operating  temperature  of  the  tether 
is  basically  independent  of  the  tether  diameter.  If  no  current  is  passed 
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through  a  tether  system  characterized  by  equal  parameters  of  absorptance 
and  emittance,  the  magnitude  of  the  temperature  fluctuations  will  thus 
I  remain  essentially  constant  as  the  diameter  increases.  This  can  be  seen 

by  rewriting  the  equation  for  the  thermal  balance  in  an  abbreviated  steady 
state  form,  as 


where  the  solar  constant  is  represented  as  S. 

The  temperature  fluctuations  due  to  ohmic  heat  generation  may  be 
evaluated  by  looking  at  the  same  equation  and  adding  the  term  for  Joule 
heating.  If  the  open  circuit  equation  listed  above  is  defined  as  a 
constant  K,  then 
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When  a  constant  current  fs  allowed  to  pass  through  tethers  of  increasing 
diameters,  the  average  temperature  would  be  expected  to  decrease  slowly 
because  the  Joule  heating  decreases  as  the  diameter  cubed. 


If  the  tether  is  required  to  generate  a  constant  power  to  the  load  for 
all  diameters,  then  the  current  level  must  be  expressed  as 
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Consequently,  the  analogous  expression  for  the  temperature  becomes 
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and  the  open  circuit  voltage  Vq^  has  been  approximated  as  v  B  L.  Thus, 
for  tethers  of  increasing  diameters  required  to  maintain  a  constant  level 
of  power  to  a  load,  the  average  temperature  will  again  begin  to  decrease 
as  the  diameter  increases. 

Figures  18A  and  18B  show  these  two  phenomena  graphically  by  portraying 
the  variation  in  average  operating  temperature  versus  increasing  tether 
diameter  for  the  cases  of  constant  current  and  current  required  to 
maintain  constant  load  power,  respectively.  As  expected,  the  average 
temperature  in  both  cases  decreases  with  increasing  diameter  when  current 
is  passed  through  the  tether. 

Tether  Length 

Similar  to  the  diameter,  the  length  of  the  tether  affects  the  amount 
of  energy  which  is  absorbed,  the  amount  which  is  radiated  away,  and  the 
magnitude  of  the  electromagnetic  force  which  is  induced  across  the  tether, 
and  hence  the  amount  of  energy  generated  as  ohmic  heat.  In  each  case,  the 
outcome  is  linearly  proportional  to  this  length.  By  referring  to  the 
equations  discussed  in  the  previous  section,  it  can  be  seen  that  a  tether 
of  arbitrary  length  would  be  expected  to  experience  basically  the  same 
changes  in  temperature  as  a  tether  of  different  length,  which  moves  in  an 
identical  orbit,  for  the  same  level  of  current.  The  overall  average 
operating  temperature  is  independent  of  variations  in  length,  as  the 
increased  amount  of  energy  absorbed  and  generated  internally  will  be 
offset  by  the  increased  amount  radiated  away.  Thus,  the  average  operating 
temperature  versus  length  for  a  tether  with  all  other  variables  remaining 
constant  is  itself  a  constant. 

Obviously,  as  the  induced  voltage  across  the  tether  grows  in  magnitude 
with  increasing  length,  more  ohmic  heat  will  be  generated  (providing  that 
the  current  level  is  not  inhibited  from  increase  and  that  all  other 
variables,  in  particular  the  load  impedance,  are  held  constant)  and  the 
operating  temperature  will  rise.  It  must  be  noted  that  the  tether  power 
increases  approximately  as  under  these  conditions,  however.  Figure  19 
illustrates  the  converse  of  this  by  comparing  the  variation  in  operating 


temperature  versus  increasing  length  for  a  tether  required  to  maintain  a 
constant  power  level  to  the  load  for  all  lengths;  when  operated  in  this 
manner,  only  progressively  decreasing  levels  of  current  are  allowed  to 
pass  through  the  system,  driving  the  average  operating  temperature  lower. 

4.1.2  Effects  of  Orbital  Parameters 
A1 ti tude 

Variation  in  the  operational  altitude  of  the  tether  system  will  affect 

both  the  orbital  velocity,  as  v  =  ^  for  a  circular  orbit,  and  the  magnetic 
field  strength,  which  is  inversely  proportional  to  the  cube  of  the  orbital 
radius.  Thus,  the  total  effect  on  the  induced  voltage  across  the  tether 
will  go  as  r'^*^.  As  the  ohmic  heat  generated  by  passing  current  is 
proportional  to  the  square  of  the  induced  voltage,  an  increase  in  altitude 
from  a  300  km  to  a  500  km  orbit  will  correspond  to  an  approximately  19% 
decrease  in  internally  generated  heat. 

Increasing  altitude  will  also  affect  the  incoming  flux  which  the 
tether  receives.  The  effect  on  solar  radiation  is  basically  negligible, 
as  the  tether  remains  in  the  sunlight  for  only  a  slightly  longer  period  of 
time,  approximately  .01%,  for  a  300  km  versus  a  500  km  orbit,  and  the 
magnitude  of  the  solar  constant  is  unchanged  for  such  a  slight  altitude 
variation.  The  other  incoming  fluxes  are  more  severely  affected,  however. 
The  albedo  flux  and  earth  radiation  flux  are  both  dependent  on  the  earth 
view  factor,  as  previously  defined  in  Section  3.2.2.  This  factor 
decreases  by  approximately  15%  from  a  300  km  to  a  500  km  orbit.  The 
aerodynamic  heating  is  dependent  on  the  atmospheric  density  and  the 
relative  velocity  at  the  operational  altitude;  the  density  decreases  by 
approximately  97%,  and  the  relative  velocity  decreases  by  approximately 
25- ,  from  a  300  km  to  a  500  km  orbit.  Thus,  the  overall  affect  on 
aerodynamic  heating  is  a  decrease  of  approximately  97%.  Radiation  from 
the  tether,  on  the  other  hand,  is  independent  of  the  altitude  of  the  orbit 
(assuming  that  the  background  temperature  at  both  300  km  and  500  km  is  the 
same) . 
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The  net  effect  then  is  that  the  tether  will  operate  at  lower 
temperatures  as  the  orbital  altitude  is  increased.  This  is  graphically 
portrayed  in  Figure  20,  showing  the  average  operating  temperature  versus 
altitude.  While  it  was  stated  previously  that  it  is  desirable  to  cause 
the  tether  to  operate  at  low  temperatures  in  order  to  minimize  ohmic 
losses  in  the  system,  this  is  correct  only  for  a  constant  altitude.  It 
must  be  remembered  that  the  magnitude  of  the  power  which  the  tether  is 
capable  of  generating  is  also  inversely  proportional  to  the  orbital  radius 
raised  to  the  seventh  power. 

Orbital  Inclination 

The  most  significant  effect  of  the  orbital  inclination  is  on  the 
magnitude  of  the  induced  voltage  across  the  tether.  At  0®  inclination, 
the  velocity  vector  of  the  tether  is  perpendicular  to  the  magnetic  field, 
thus  maximizing  the  induced  voltage.  Similarly,  in  a  polar  orbit,  or  an 
inclination  of  90®,  the  velocity  is  parallel  to  the  magnetic  field  and  no 
voltage  is  induced.  Because  the  induced  voltage  is  calculated  as  the 
cross  product  of  y  and  B,  the  magnitude  of  the  induced  voltage  varies  as 
cos  i  (assuming  a  rigid  tether  aligned  with  the  vertical). 

The  operating  temperature  of  the  tether  is  relatively  insensitive  to 
changes  in  the  orbital  inclination,  with  none  of  the  energy  inputs  being 
significantly  affected.  With  regard  to  solar  radiation,  the  orbital 
inclination  primarily  changes  the  amount  of  surface  area  of  the  tether 
which  is  exposed  to  the  incoming  flux,  and  hence  how  much  energy  is 
instantaneously  absorbed;  however,  the  average  over  the  course  of  an  orbit 
is  independent  of  the  inclination.  With  regard  to  earth  radiation  and 
albedo,  it  primarily  changes  the  amount  of  incoming  energy  because  of 
variations  in  cloud  cover,  surface  conditions  of  snow  and  ice,  and  water 
covered  areas.  In  low  inclination  orbits,  which  are  the  only  feasible 
ones  for  a  tether  system  because  of  the  decrease  in  induced  voltage  at 
high  inclinations,  typical  variations  between  0®  and  28.5®  are 
approximately  a  6%  increase  for  earth  radiation  and  albedo  flux,  and  a  3% 
increase  for  aerodynamic  heating. 
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Longitude  of  the  Line  of  Nodes 


It  would  be  expected  that  the  longitude  of  the  line  of  nodes  of  the 
orbit  would  have  little  effect  on  the  temperature  of  the  tether.  The  only 
variable  which  it  affects  is  the  flux  due  to  solar  radiation  by  varying 
the  magnitude  of  the  surface  area  of  the  tether  which  is  exposed  to  the 
incoming  flux.  Typically,  shifting  the  line  of  nodes  by  90®,  from  a  line 
perpendicular  to  the  earth-sun  radius  to  a  line  parallel  to  this  radius, 
will  cause  a  decrease  of  approximately  68%  in  the  solar  radiation 
absorbed. 

It  can  be  seen  from  Figure  22  that  even  this  decrease  in  the  solar 
flux  has  little  effect  on  the  operating  temperature,  acting  only  to 
slightly  decrease  the  average  temperature  and  limiting  the  amplitude  of 
modulation.  Both  cases  shown  here  are  with  the  ecliptic  and  equatorial 
planes  aligned.  Figure  21A  represents  the  tether  in  a  circular  orbit  with 
the  line  of  nodes  oriented  perpendicular  to  the  earth-sun  line;  it  is  this 
orientation  which  exposes  the  tether  to  the  largest  incoming  flux  during 
the  course  of  its  orbit.  The  average  temperature  in  this  orientation  is 
approximately  215°  K,  with  temperature  ranges  over  35®.  Figure  21B 
represents  an  orbit  with  the  line  of  nodes  shifted  by  90®  so  that  it  is 
aligned  to  the  earth-sun  line.  Now  the  average  temperature  is  slightly 
lower,  as  expected,  at  approximately  205®  K  with  a  temperature  range  of 
20®. 

Seasonal  Variations 


Seasonal  variations,  meaning  fluctuations  due  to  the  position  of  the 
earth  in  its  orbit  about  the  sun,  are  generally  of  two  types.  First,  the 
magnitude  of  the  solar  flux,  as  measured  at  the  earth,  varies  because  of 
the  eccentricity  of  the  earth's  orbit.  Second,  the  tilt  of  the  earth's 
rotational  axis  to  the  ecliptic  causes  the  angle  between  the  earth-sun 
line  and  the  equator  to  vary  between  ±  23.5®.  The  first  of  these  affects 
all  of  the  external  fluxes  except  for  aerodynamic  heating,  while  the 
second  affects  only  the  solar  radiation  flux.  However,  the  magnitude  of 
the  variations  is  small  enough  so  that  it  would  be  expected  that  the 
operating  temperature  of  the  tether  will  not  be  altered  significantly. 
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Typically,  the  solar  flux  varies  by  approximately  3%  throughout  the  year, 
while  the  sun  angle  causes  variations  in  the  solar  radiation  flux  of 
approximately  9%. 

Since  the  variations  in  the  solar  flux  are  much  less  than  they  were 
for  changes  in  the  line  of  nodes,  it  is  anticipated  that  the  operating 
temperature  of  the  tether  would  be  affected  very  little  by  seasonal 
changes.  Figures  22  show  this  to  be  the  case.  The  three  figures 
presented  represent  the  tether  at  the  summer  solstice  (approximately 
aphelion),  at  the  winter  solstice  (approximately  perihelion),  and  at  the 
vernal  equinox.  Only  very  small  variations  in  the  average  temperature  and 
the  range  of  temperatures  are  discernible  between  the  three  graphs. 

4.2  VARIATIONS  OF  ORBITAL  ELEMENTS 

The  manner  in  which  the  tether  is  operated,  i.e.  as  a  power  generator, 
thruster,  or  a  combination  of  the  two,  may  affect  the  elements  of  the 
orbit  in  which  it  begins  its  operation.  This  is  because  the  production  of 
the  induced  electromagnetic  drag  and/or  thrust  force  may  not  always  be 
constant  throughout  the  orbit,  nor  may  it  always  lie  in  the  orbital  plane. 
Indeed,  it  is  typical  that  neither  of  these  two  conditions  will  be  true. 

This  examination  will  use  the  method  described  in  Section  3.3  to 
analyze  the  variation  of  the  orbital  elements  for  the  cases  of  the  tether 
orbiting  but  not  passing  any  current,  for  the  pure  generator  mode,  for  the 
pure  thruster  mode,  and  finally  for  the  case  of  the  tether  generating 
power  when  in  the  earth's  shadow  and  thrusting  when  in  sunlight. 

The  orbital  elements  which  are  looked  at  are  the  semi -major  axis, 
which  corresponds  to  the  energy  in  the  orbit,  the  eccentricity,  the 
longitude  of  the  line  of  nodes,  the  orbital  inclination,  and  the  argument 
of  pericenter.  These  are  calculated  as  described  in  Section  3.3.2.  Thus, 
the  variation  of  the  semi-major  axis  and  the  eccentricity  are  affected 
only  by  the  components  of  the  disturbing  acceleration  in  the  orbital 
plane,  the  variation  of  the  longitude  of  the  line  of  nodes  and  the  orbital 


inclination  are  affected  only  by  the  out-of-plane  component,  while  the 
argument  of  pericenter  varies  with  both  in  and  out-of-plane  forces.  These 
may  be  conveniently  illustrated  by  writing  these  variational  equations  as 


^  _  r  sin  d  ■ 
dt  ”  h  sin  i  -h  *  -d 


^  =  r  cos  6  . 
dt  h  -h  ‘  -d 


cos  i 


dt 


For  this  analysis,  initial  conditions  were  arbitrarily  selected  for  a 
400  km  orbit  of  28.6®  inclination,  with  the  line  of  nodes  perpendicular 
to,  and  the  line  of  apsides  parallel  to,  the  earth-sun  line.  The  orbit 
was  specified  as  near-circular,  with  a  finite  eccentricity  of  .01.  The 
tether  itself  was  specified  as  an  aluminum  conductor,  3  mm  in  diameter  and 
20  km  long.  The  total  mass  of  the  system,  tether  and  sub-satellite,  was 
specified  as  25000  kg. 


4.2.1  Open  Circuit  Operation 

In  the  open  circuit  operation,  no  current  is  allowed  to  pass  through 
the  tether;  hence,  no  electromagnetic  force  is  produced.  Therefore,  the 
disturbing  acceleration  acting  on  the  tether  is  composed  only  of 
aerodynamic,  gravitational,  libration,  and  solar  radiation  specific 
forces,  and  the  tether  system  is  effectively  just  another  earth-orbiting 
satellite.  Each  element  of  the  disturbing  acceleration  has  both  in  and 
out-of-plane  components,  with  the  exception  of  the  libration  specific 
forces  which  act  only  in-plane  because  of  our  first  order  approximation. 
As  described  in  Section  3.3.1,  the  dominant  element  of  these  is  the 
contribution  due  to  the  earth's  oblateness.  The  other  components  cause 
intra-orbital  fluctuations  in  the  orbital  elements,  and  can  have  definite 
long-term  effects  on  the  orbit,  but  in  the  short-term  generally  do  not 
appreciably  perturb  the  shape  of  the  orbit. 
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To  make  an  initial  estimate  of  the  effect  of  this  disturbing 
acceleration  on  the  orbital  elements,  it  is  convenient  to  use  Gauss' 
planetary  equations  (the  derivation  of  which  is  included  in  Appendix  C, 
for  a  coordinate  frame  with  radial,  tangential,  and  normal  components)  and 
linearize  them  by  assuming  that  the  eccentricity  of  the  orbit  is  small. 
When  this  done,  Gauss'  variational  equations  become 
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The  representation  of  the  disturbing  acceleration  due  to  the  earth's 
oblateness  is  indicated  in  Appendix  B.  The  vector  components  for  this 
axially  symmetric  representation  may  be  written  as 
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When  these  are  substituted  into  the  variational  equations,  the  mean 
variational  rates  of  the  orbital  elements  may  be  calculated.  Doing  this, 
with  values  for  the  orbit  of  our  system,  indicates  that  the  inclination, 
semi -major  axis,  and  eccentricity  remain  essentially  constant  throughout 
an  orbit,  while  the  longitude  of  the  line  of  nodes  and  the  argument  of 
peri  center  vary  as 
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<  ^  y  =  -7.1  degrees/day 
and 

<  ^  >  =  11.5  degrees/day 

These  average  secular  variations  in  Q  and  u  are  classically  referred  to  as 
the  regression  of  the  node,  and  the  precession  of  the  argument  of  perigee, 
respectively.  For  an  earth  orbiting  satellite,  the  magnitudes  of  these 
average  variations  may  also  be  calculated  directly  from  the  disturbing 
function  associated  with  the  earth's  gravitational  field  and  expressed  as 
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which  corroborate  the  estimates  calculated  previously. 


Figures  23  depict  the  variations  of  the  orbital  elements  over  an  8  day 
period.  As  can  be  seen,  the  mean  variations  are  essentially  as 
anticipated  above,  with  intra-orbital  fluctuations  superimposed  onto  the 
mean.  The  slight  decrease  in  the  eccentricity  and  in  the  semi-major  axis 
with  time  may  be  attributed  to  the  effects  of  aerodynamic  drag  as  it  takes 
energy  from  and  begins  to  circularize  the  orbit.  All  satellites  in  orbits 
identical  to  the  tether  system  will  experience  variations  in  their  orbit 
exactly  as  the  tether  in  the  absence  of  the  electromagnetic  force.  When 
current  is  permitted  to  flow  through  the  tether,  the  variations  caused  by 
it  will  be  superimposed  on  the  base  variations. 


4.2.2  Pure  Mode  Operation 

As  stated  above,  the  only  difference  between  pure  mode  operation  and 
open  circuit  operation  is  the  addition  of  the  induced  electromagnetic 
force  to  the  disturbing  acceleration.  When  operated  in  only  a  single 
mode,  the  components  of  the  electromagnetic  force  are  cyclic  with  the  true 
anomaly.  Assuming  that  the  tether  remains  vertical,  these  may  be 
approximated  as 
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where 


=  0 

at  =  B  (C  +  cos  2f) 
ap  =  D  sin  f 

B  =  +3.635  X  10“6  ^2  (positive  for  thrusting) 
C  =  15.504 
D  =  1.735  X  10-6 


It  can  be  seen  that  the  electromagnetic  force  dominates  in  the  tangential 
direction,  while  the  oblateness  effects  continue  to  dominate  in  the  radial 
and  normal  directions.  When  these  accelerations  are  added  to  those  which 
are  present  before  the  current  is  passed,  and  the  sum  is  substituted  into 
the  variational  equations  as  before,  the  mean  variational  rates  which  are 
to  be  expected  when  the  tether  is  operated  in  a  single  mode  are 


<  ^  >  =  +  .01  degrees/day  (negative  for  thrusting) 

<  ^  >  =  -7.1  degrees/day 

^  Ht  ^  ^  km/day  (positive  for  thrusting) 

<3f  >  =  ° 


<  ^  >  =  11.5  degrees/day 

Not  surprisingly,  the  elements  all  essentially  vary  as  they  do  for  an 
ordinary  satellite  with  the  exception  of  the  semi-major  axis,  which 
increases  as  energy  is  added  when  the  tether  is  operated  as  a  thruster  and 
decreases  when  it  is  operated  as  a  generator,  and  the  orbital  inclination, 
which  varies  in  the  opposite  manner.  The  eccentricity  also  actually 
varies  slightly,  but  not  to  noticeable  magnitudes. 


Figures  24  and  25  verify  these  expectations.  Reviewing  the  graphs  of 
the  variations  of  the  orbital  elements  for  a  pure  generator  operation  as 
depicted  in  Figures  24,  and  a  pure  thrusting  operation  as  depicted  in 
Figures  25,  it  can  be  seen  that  the  orbital  inclination  does  display  a 
finite  change  with  time,  and  the  semi-major  axis  and  perigee  altitude 


58 


■■T.vv 


- 


change  dramatically  in  response  to  the  operational  mode.  As  before,  there 
are  intra-orbital  fluctuations  superimposed  onto  the  mean  variations;  the 
magnitude  of  the  amplitudes  of  these  fluctuations  does  not  appear 
significantly  different  from  those  calculated  before  the  emf  was  added. 
No  significant  variations,  beyond  those  which  were  anticipated,  are 
apparent  in  the  other  orbital  elements  over  the  8  day  period. 

4.2.3  Mixed  Mode  Operation 

Mixed  mode  operation  is  characterized  by  alternately  causing  the 
tether  to  generate  and  thrust  throughout  the  course  of  its  orbit,  the 
purpose  of  which  is  to  maintain  the  average  energy  of  the  orbit,  and  thus 
the  semi-major  axis,  at  an  approximately  constant  level.  The  result  of 
this  is  a  continual  variation  in  the  direction  of  the  disturbing 
acceleration.  This  acts  to  intensify  the  variations  in  the  orbital 
elements  as  the  tether  completes  an  orbit. 

Since  the  disturbing  acceleration  due  to  the  induced  electromagnetic 
force  varies  throughout  the  orbit,  and  not  as  a  function  of  the  true 
anomaly,  an  analytical  approach  via  Gauss’  equations  as  before  does  not 
seem  appropriate.  Although  the  induced  forces  may  be  generally  considered 
as  being  continually  applied  throughout  the  orbit,  it  is  convenient  for  a 
simple  analysis  to  consider  them  as  impulsive  forces  applied  at  the 
perigee  and  apogee  of  the  orbit.  Thus,  if  the  operating  mode  being 
considered  is  one  in  which  the  energy  of  the  orbit  is  being  held  constant 
by  generating  power,  or  extracting  energy  from  the  orbit,  in  the  earth's 
shadow,  and  injecting  an  equal  amount  of  energy  into  the  orbit  by 
thrusting  when  in  sunlight,  then  we  can  consider  this  as  a  finite  force 
applied  in  the  direction  of  velocity  on  the  sunlight  side  of  the  orbit, 
and  an  equal  but  opposite  force  applied  on  the  shadow  side.  If  the  orbit 
is  initially  specified  as  being  characterized  by  a  small,  finite 
eccentricity  vector  pointing  toward  the  sun  and  parallel  to  the  earth-sun 
line,  then  this  operating  mode  is  similar  to  that  of  applying  an  impulsive 
A  V  in  the  direction  of  the  orbital  velocity  at  perigee,  and  applying  a 
second,  but  equal,  one  opposite  the  orbital  velocity  at  apogee.  The 
maneuver  at  perigee  would  be  expected  to  add  energy  to  the  orbit,  and 
consequently  increase  the  semi-major  axis  and  apogee  height,  while  the 
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apogee  maneuver  would  subtract  energy,  and  consequently  decrease  the  semi- 
major  axis  and  the  perigee  height.  The  cumulative  effect  of  this 
alternately  thrusting  and  generating  is  to  increase  the  eccentricity  of 
the  orbit  as  the  altitude  of  perigee  decreases  and  the  altitude  of  apogee 
increases.  The  other  orbital  elements  are  unaffected  by  this  type  of 
operation,  except  as  they  are  influenced  by  the  remainder  of  the 
components  of  the  disturbing  acceleration,  just  as  was  described  in  the 
previous  two  sections. 

Figures  26  graphically  illustrate  this  phenomenon  by  portraying  the 
variations  of  each  of  the  orbital  elements  over  an  8  day  period  when 
operated  in  the  mixed  mode  described  above.  As  can  be  seen  the  variations 
in  the  argument  of  perigee  and  in  the  longitude  of  the  line  of  nodes 
remain  essentially  as  was  determined  for  the  open  circuit  mode  of 
operation.  The  variation  in  the  semi -major  axis  may  be  attributed  to  an 
imbalance  between  the  amount  of  time  that  the  tether  is  caused  to  generate 
power  versus  the  amount  that  is  caused  to  thrust.  Generally,  it  may  be 
observed  that  it  experiences  no  significant  change  since  the  energy  of  the 
orbit  remains  essentially  constant.  The  eccentricity,  and  unexpectedly 
the  orbital  inclination,  however,  increase  relatively  rapidly  over  the 
time  period  considered.  If  the  rate  of  eccentricity  change  is 
extrapolated  form  this  8  day  period,  it  indicates  that  after  approximately 
85  days  the  tether  will  descend  to  a  level  at  which  the  perigee  will  be 
low  enough  to  cause  the  tether  to  re-enter  the  atmosphere.  Fortunately, 
this  does  not  appear  to  be  precisely  the  case,  as  will  be  shown  shortly; 
however,  this  type  of  operation  will  still  be  prohibitive  in  and  of 
i tsel f . 

Using  the  type  of  simple  analytical  approach  as  before  to  develop  a 
solution  to  this  problem  erroneously  leads  one  to  a  mode  of  operation  in 
which  the  tether,  initially  inserted  into  an  orbit  which  has  its  perigee 
opposite  the  sun,  still  generates  power  throughout  the  shadow  period,  but 
thrusts  to  make  up  the  orbital  energy  loss  only  in  the  vicinity  of  the 
line  of  nodes.  While  it  would  appear  that  the  effect  of  this  operating 
technique  is  to  continually  lower  the  apogee  height  while  maintaining  the 


semi -major  axis  constant  over  the  course  of  each  orbit,  when  this 
technique  is  simulated  the  eccentricity  of  the  orbit  still  displays  a 
rapid  increase,  as  before,  and  the  system  rapidly  gets  low  enough  to  re¬ 
enter  the  atmosphere.  In  fact,  this  increase  in  eccentricity  was  evident 
in  every  other  scheme  which  could  be  devised  for  operating  the  tether  in 
mixed  mode.  This  is  evidently  a  problem  for  which  there  is  no  intrinsic 
solution;  if  the  tether  is  to  be  operated  in  this  manner,  compensation  for 
the  eccentricity  increase  will  have  to  be  externally  provided. 

If  the  power  to  mass  ratio  of  the  system  is  increased,  by  increasing 
the  current  level  and/or  decreasing  the  system  mass,  the  saturating 
effects  of  the  eccentricity  variation  may  be  observed.  When  this  is  done 
with  the  system  current  and  mass  specified  at  different  levels,  it  can  be 
seen  in  Figure  27  that  the  increase  in  the  eccentricity  is  actually  more 
sinusoidal  than  linear,  with  a  maximum  amplitude  of  .08  and  a  period  of 
approximately  165  days  when  referred  back  to  the  original  system 
specifications.  The  cause  of  this  saturating  effect  is  due  to  the 
simultaneous  variations  in  the  longitude  of  the  line  of  nodes  and  the 
argument  of  perigee,  caused  by  the  earth's  oblateness.  However,  as  an 
eccentricity  in  excess  of  approximately  .04  insures  that  the  system  will 
re-enter  the  atmosphere  before  the  orbit  re-ci rcularizes,  the  magnitude  of 
this  amplitude  still  makes  the  use  of  the  mixed  mode  of  operation 
prohibitive  for  a  tether  system  which  utilizes  a  constant  current  without 
an  external  system  dedicated  to  compensate  for  the  eccentricity  increase. 


5.  CONCLUSION 


5.1  SUMMARY 

Within  the  constraints  of  the  assumptions  which  were  made  to  conduct 
this  analysis,  several  conclusions  may  be  drawn.  Some  of  them  are 
perfectly  obvious,  while  others  are  somewhat  unexpected. 

First,  the  tether  may  be  operated  in  an  uninsulated  fashion  as  well  as 
an  insulated  one.  The  insulated  mode  must  face  the  problems  of  arcing  and 
dielectric  degradation  due  to  the  tether's  operational  environment,  in 
addition  to  leakage  losses  of  comparatively  large  magnitudes  through  any 
holes  which  are  present  in  the  insulating  jacket.  The  uninsulated  mode 
must  face  current  leakages  along  the  tether's  length  on  the  order  of  10% 
plus  some  minor  degradation  of  the  conductor  itself  due  to  sputtering. 
The  surprising  result  is  that  the  tether  need  not  actually  be  fully 
insulated  against  the  surrounding  plasma  in  order  to  operate;  however, 
this  type  of  operation  would  be  done  at  the  expense  of  the  current  leakage 
as  described  in  Section  3.1,  and  over  wider  operating  temperature  ranges 
oecause  of  the  relatively  shiny  conductor  surface  as  described  in  Section 
4.1.1.  The  use  of  anodized  aluminum  as  the  dielectric  coating  for  the 
conducting  material  of  the  tether  appears  to  be  an  attractive  solution  to 
both  types  of  problems.  It  is  characterized  by  high  dielectric  strength, 
durability,  high  emittance  and  low  absorptance.  Its  use  as  an  insulator 
certainly  merits  further  theoretical  and  experimental  analysis. 

The  operating  temperature  of  the  tether  is  determined  primarily  by  its 
values  of  absorptance  and  emittance.  It  is  not  significantly  affected  by 
tether  length,  nor  by  tether  diameter  once  a  certain  minimum  is  reached. 
To  operate  at  low  temperatures,  and  hence  at  low  tether  resistances  and 
higher  efficiencies,  higher  ratios  of  emittance  to  absorptance  must  be 
achieved  by  properly  selecting  the  surface  material  of  the  tether.  In 
general,  the  operating  temperatures  range  by  ±  20  degrees  about  the 
average,  determined  by  the  surface  characteristics,  over  the  course  of 
each  orbit.  Changes  in  the  orbital  elements,  for  near  circular  orbits, 
and  seasonal  variations,  as  the  earth  orbits  the  sun,  do  not  generally 


change  the  average  operating  temperature,  but  only  modulate  the  intra¬ 
orbital  temperature  range.  As  would  be  expected,  tether  resistance 
basically  follows  temperature  variations,  while  current  flow  and  induced 
force  basically  follow  the  magnitude  of  the  induced  potential  across  the 
tether.  Efficiency  is  affected  by  both  temperature  and  current,  generally 
being  inversely  proportional  to  current. 

The  type  of  operation  in  which  the  tether  system  is  employed  can 
affect  the  variation  of  the  orbital  elements  over  time.  The  open  circuit 
mode,  generator  mode,  and  thruster  mode  present  no  surprises.  As 
anticipated,  the  only  variations  present  when  no  current  flow  is  permitted 
are  those  due  primarily  to  the  earth's  oblateness,  causing  precession  of 
the  argument  of  pericenter  and  regression  of  the  line  of  nodes,  but  only 
negligibly  small  effects  on  the  other  orbital  elements.  Also  as 
anticipated,  the  generator  mode,  in  addition  to  the  variations  present  in 
the  open  circuit  mode,  causes  the  semi-major  axis  to  dramatically  decrease 
as  the  orbital  energy  of  the  tether  is  transformed  into  on-board  power. 
Similarly,  the  semi -major  axis  increases  rapidly  in  the  thruster  mode  as 
energy  is  added  to  the  orbit  by  the  tether's  operation.  Surprisingly,  in 
both  modes  the  orbital  inclination  make  small  but  detectable  changes, 
increasing  in  the  generator  mode  and  decreasing  when  thrusting.  Beyond 
this,  the  other  orbital  elements  show  no  other  significant  variations 
beyond  those  evident  in  the  open  circuit  mode.  In  mixed  mode,  however, 
while  the  semi-major  axis,  and  hence  the  orbital  energy,  can  be  held 
constant  when  averaged  over  an  orbit,  it  appears  impossible  when  operating 
with  a  constant  current  to  preclude  an  increase  in  orbital  eccentricity  to 
the  extent  that  perigee  is  soon  low  enough  for  the  tether  to  re-enter  the 
atmosphere.  This  effect  may  be  mitigated  by  operating  with  a  low  system 
power-to-mass  ratio,  which  slows  the  rate  of  this  increase  to  what  are 
probably  tolerable  levels. 

5.2  RECOMMENDATIONS 

All  of  this  analysis  has  been  on  the  basis  of  several  simplifying 
assumptions.  Essentially,  two  types  were  made:  those  assumed  because  it 
is  known  that  the  assumption  is  a  good  one,  either  because  parameters 
neglected  are  small,  linearization  was  desired  for  simplicity,  or  models 
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used  have  been  shown  to  be  accurate  representations;  and  those  assumed 
because  to  date  there  is  no  definitive  information  available  which  can  be 
used  to  produce  a  more  accurate  model  for  the  analysis.  The  former  group 
includes,  among  others,  the  assumptions  of  the  tether  modelled  as  a  rigid 
body  composed  of  point  masses,  and  the  mathematical  models  of  the  earth's 
magnetic  field,  mass  distribution,  ionospheric  composition,  and  air 
density.  The  latter  group  primarily  includes  those  made  to  produce  the 
simplified  electrical  circuit  to  model  the  tether  system,  and  to  neglect 
any  dynamic  effects  of  the  tether's  operation,  such  as  oscillations  and 
deployment.  Because  of  the  lack  of  any  actual  experimental  or  flight 
data,  most  of  what  is  available  in  these  areas  is  theory;  its  accuracy  is 
the  subject  of  significant  debate.  It  is  on  these  subjects  that  more 
extensive  research,  coupled  with  actual  flight  testing,  should  be 
dedicated  to  determine  if  the  model  which  as  been  used  here  is  adequate  or 
must  be  modified  to  produce  accurate  results. 

The  main  assumptions  which  have  been  made  regarding  the  system's 
circuit  path  are  the  following:  the  impedance  of  the  ionosphere  is 
negligible;  the  potential  drop  across  each  contactor  which  is  required  to 
permit  sufficient  current  to  flow  through  the  circuit  is  relatively  small, 
on  the  order  of  100  volts  as  assumed  here,  compared  to  the  total  induced 
potential  across  the  tether;  and  the  properties  of  anodized  aluminum  are 
such  that  the  tether  would  be  essentially  impervious  to  sputtering  and 
current  leakage,  while  possessing  the  required  characteristics  of 
absorptance  and  emittance.  Those  made  regarding  dynamics  are  that  the 
deployment  of  the  tether  from  the  base  satellite  may  be  accomplished  with 
little  difficulty,  and  the  tether  will  remain  deployed  permanently;  and 
that  the  dynamics  of  the  tether  motion,  particularly  resonance  with  the 
out-of-plane  oscillation  frequency,  do  not  pose  a  problem  to  the  system's 
operation.  If  the  impedance  of  the  ionosphere  is  much  more  than  the  few 
ohms  which  were  assumed  in  the  discussion  of  Section  2.2.1,  then  the 
efficiency  of  the  system,  and  hence  its  ability  to  compete  with  other 
power  generation  systems,  will  drop  below  what  may  be  considered  an 
acceptable  level.  Similarly,  if  the  potential  drops  across  the  anode 
and/or  cathode,  which  must  be  developed  in  order  to  pass  sufficient 


current  through  the  tether  to  reach  the  desired  levels  of  power 
generation,  become  so  large  that  they  also  drive  the  overall  efficiency  of 
the  system  below  acceptable  levels,  then  the  system  will  not  be 
competitive.  Obviously,  if  the  tether  itself  cannot  be  deployed  because 
of  dynamical  difficulties,  or  cannot  be  kept  in  operation  without 
excessive  oscillations,  then  the  system  will  be  physically  impossible  to 
employ.  Finally,  if  a  suitable  dielectric  cannot  be  found  which  possesses 
the  necessary  qualities  of  high  emittance,  low  absorptance,  high 
dielectric  strength,  and  resistance  to  the  operational  environment,  then 
it  will  be  necessary  to  accept  the  additional  losses  incumbent  on 
operating  as  an  uninsulated  conductor. 

Manifesting  an  experiment  on  a  shuttle  flight,  as  is  currently 
scheduled  in  the  form  of  the  TSS,  should  provide  answers  to  most  of  these 
questions.  On  smaller  scales,  the  magnitude  of  the  impedance  of  the 
ionosphere  and  contactors  could  be  determined,  as  well  as  the  motions  of 
the  tether  once  it  is  deployed  and  in  operation.  Deployment  schemes  would 
have  to  be  exercised,  as  well  as  the  additional,  and  probably  greater, 
problems  involved  in  retrieval.  Prior  to  actual  flight  testing,  further 
research,  both  experimental  and  theoretical,  should  be  devoted  to  the 
feasibility  of  using  anodized  aluminum  as  a  dielectric  coating  for  the 
conductor,  and  to  the  development  and  testing  of  newer  and  more  efficient 
contactor  configurations. 
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APPENDIX  A 


Computer  Calculation  of  Thermal  Variations 

The  thermal  balance  of  the  tether  system,  and  the  variations  in  the 
temperature,  resistance,  current  level,  and  efficiency  of  the  system,  were 
evaluated  with  the  mathematical  model  as  discussed  in  Section  3.2.  This 
appendix  includes  the  computer  program  which  implements  this  model.  It 
employs  an  Euler  forward  code  for  the  numerical  integration  of  the 
temperature  equation  with  a  time  step  of  60  seconds.  Although  this  is  a 
relatively  inaccurate  integration  method,  the  purpose  of  the  analysis 
which  this  program  implements  is  to  show  general  trends  in  the  variation 
of  system  values.  Any  numerical  errors  which  develop  over  the  course  of 
the  integration  may  be  confidently  ignored. 

The  values  of  the  system  parameters  used  in  the  program  are  tabulated 
below. 


Physical  Constants 

Stefan-Bol tzman  Constant, 

Solar  Constant,  S 
Atmospheric  Density  at  400  km, 

Earth  Radiation 
Earth  Mean  Radius,  r^ 

Earth  Gravitational  Constant,  p 


5.669  X  10-8  w/m2/®K4 
1368  W/m2 

1.49  X  10“^^  kg/m8 
237  W/m2 
6378.1  km 

3.986  X  105  km3/s2 


Material  Properties  of  Aluminum 


Specific  Heat,  C 

Temperature  Coefficient  of  Resistance, 
Temperature  Coefficient  of  Expansion, 
Temperature  Coefficient  of  Density,  ®ni 
Resistivity  at  293.15®K 
Density  at  273.15“K 


878.4  J/kg/“K 
.00424  fi/®K 
12.9  m/m/®K 

6.9118  X  10-5  kg/m3/®k 
2.857  X  10-8  n 
2710  kg/m3 


SH 


esian  Parameters  of  Tether  &  Orbit 


Tether  Length,  L 

20  km 

Velocity,  v 

7669  m/s 

Longitude  of  Line  of  Nodes 

0® 

Sun  Zenith  Angle 

0® 

Orbital  Inclination 

28.5® 

Earth  View  Factor,  F 

.31413 

Absorptivity 

.15 

Emissivity 

.85 

Earth  Albedo  Factor 

.28 

Contactor  Potential  Drop 

100  V 

Orbital  Altitude 

300  km 

The  program  presumes  the  tether  to  be 
prescribed  orbital  altitude. 

in  a  circular  orbit  at  the 
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C  THIS  PROGRAM  COMPUTES  THE  VARIATION  IN  TEMPERATURE  OF  A 
C  TETHER  WHICH  IS  OPERATED  IN  MIXED  MODE.  THE 

C  THRUSTING  PHASE  MAY  BE  DONE  WHILE  MAINTAINING 

C  EITHER  CURRENT,  VOLTAGE,  OR  FORCE  CONSTANT. 

C 

C  PARAMETERS  USED  ARE  FOR  A  300  KM  ORBIT. 

C 

C  DATA  FILES  ARE  NUMBERED  SEQUENTIALLY  FOR  TEMPERATURE, 

C  RESISTANCE,  CURRENT,  FORCE  8.  EFFICIENCY. 

C 

C 

c 

REAL  I TER . I NC , L , MASS , LONG , L I TE , LENGTH 
DIMENSION  T( 2000 ) , VB<2000 ) , VMAG( 2000 ) . DOT( 2000 ) 

C 

C 

C  SPECIFY  SYSTEM  PARAMETERS 

C 

c 

PARAMETER( P I =3 . 1 4 1 5927 , RH0R=27 10., AEXP= 1 2 . 9D-6 , LONG=0 . ) 
PARAMETERC  ALFA= , 00424 , F= . 31 41 3 , B0LT=5 . 669D-8 , RATE= 1 . 1 569D-3 ) 
PARAMETER(C=878.4,RESTR=2.857D-8,DELT=60 . ,ERR=1 .D5,BETA=0 . ) 
PARAMETERC  V=7725 .8,1 TER= 1 727 . , I NC=28 . 5 , CRH0=6 . 9 1 1 8D- 5 ) 
PARAMETER( VCONT=200 . , R0RB=6678 . 1 , A= . 1 5 , EM= . 85 , SOL= 1368.) 
PARAMETER(ALB=.28,ERAD=237. ,RE=6378. 1 , LENGTH=2 . D4 ) 

C 

>C 

C  READ  INDUCED  VOLTAGE  FILE 

C 

C  (DATA  FILE  READ  HERE  IS  THE  ONE  PRODUCED  BY  SAO 

C  USING  NEWMAG  CODE) 

C 


CO  J= 


?EACc-iC.».£ND=50)T(  J)  .VB;  J;  .VMAG;.,:: 

vBc.:  )=VB(  J)-voo';t 


.  I  .  n  r  '  1 
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C 


WRITEC6.#)  IS  ITERATION  FOR  ENEPGV  BALANCE  TO  BE 
INCLUDET-? 


THIS  ITERATION  SERVES  TO  DETERMINE  THE  APPROPRIATE 
DIAMETER  TETHER  TO  MAINTAIN  A.  DESIRED  ^EVEl  OF  CD'PPENT. 
VOLTAGE  OR  FORCE  DEPENDING  ON  THRUSTER  OPERATING  MODE 


WRITER  6,  <► )  (YES  =  .;  NO  =  2: 
WPITLce.^t)  CHOICE  = 

READ  '(ID  .KKK 


ENTER  SPECIFICS  OF  TETHER  SY’STEM  S  PHYSICAL  SIZE.  POWER, 
LEVEL.  AND  INITIAL  POSITION  IN  ORBIT 

WRITE<6.#)' ENTER  TETHER  DIAM  (MM)  AND  PHASE  ANGLE  (RAD)' 
WRITE(6,«)'DIAMETER  =  ' 

READ  '(F12.6)DXD 
WRrTE<6,*)'PHASE  ANGLE  = 

READ  '(F6.1)MHETA0 

WRITE<6,«)  ENTER  LOAD  POWER  OF  THE  SYSTEM  (KW)' 

READ  '(F7.3)',P 

SPECIFY  TYPE  OF  THRUSTER  OPERATION 

WRITEce,#)' IS  THIS  RUN  FOR  CONSTANT  CURRENT.  CONSTANT 
+  VOLTAGE.  OR  CONSTANT  THRUST?' 

WRITE(6,#)' (CURRENT  =  1:  VOLTAGE  =  2;  THRUST  =  3) 
WR:TE(6.«)  CHOICE  - 
READ  (II)  ,NN 


C  OPEN  DATAFILES 

I . 


10? 


V 


IF(NN.EQ.1)THEN 

0PEN(51 .FiLE='F0R05l  ,STATUS=  NEW  ) 
OPEN(52.FILE=-FOR052  ,STATUS='NEW  ) 
0PEN(53.FILE=-FOR053- .STATUS='NEW  ; 

0PEN(54.FILE-  FOR054  .STATUS=^NEW  ) 

0PEN(55.FILE=  FOR055  .STATUS=' NEW- ) 

END  IF 

IF(NN.EQ.2)THEN 

0PENC61,FILE=  FOROer  .STATUS=  NEW  ) 
OPEN(62.FiLE='FDR062' . STATUS- ' NEW '  ) 
OPEN(63.FILE=FOR063' .STATUS='NEW  ) 

0PEN(64.FILE-  FOR064, STATUS- 'NEW  ; 

0PEN(65,FILE=  F0R065- .STATUS-NEW  ) 

END  IF 

IF(NN.E0.3)THEN 

0PEN(71 .FILE-  FOR071- .STATUS-NEW' ) 

0PEN(72.FILE-  FOR072' .STATUS-  NEW' ) 
0PEN(73.FILE='F0R073' .STATUS-  NEW' > 

OPEN (74, FILE-  FORD 74' .STATUS- ' NEW' ) 
OPEN<75,FILE='FOR075' .STATUS-NEW' ) 

END  IF 

ENTER  RUN  VARIABLES 

(MAY  SPECIFY  DESIRED  LEVEL  OF  VOLTAGE.  CURRENT  OR  FORCE 
TO  BE  MAINTAINED  AS  A  CONSTANT  WHEN  TETHER  IS  THRUSTING) 

IF(NN.EQ.2)THEN 

WRITE(6,#) 'ENTER  SOLAR  ARRAY  VOLTAGE  (KW)' 
WRITE(6.»)'V0LTAGE  =  ' 

READ  '(F7.1)'.V0C 
GO  TO  80 

END  IF 

IF(NN.E0.3)THEN 

WPITE<6.*)  ENTER  CONSTANT  THRUST  (N) 

WPITE(6,*:  THPUST  = 

READ  (F8.3)  .THRUST 
GC  TC  80 

nIiD  1  F 

wr;tE(6.*)  enter  no>';nal  current  ;rs?s; 

*  r' .  i  t  <  6 ,  *  )  C R  Rr,N  i  - 


ooooooo  ooooooo 


«««««##«*# 

SPECIFY  INITIAL  CONDITIONS 

»»##«»#####«»#»*»«***##»«###»#»######### 

D=XD/1000. 

XINC=INC*PI/180. 

ZEN=BETA*PI/180. 

PREV=C . 

SHADOW=ACOS<SORT( 1 .-(RE/R0RB)»#2. )/C0S(ZEN) ) 

THETA=THETA0 
TEMP=205. 

TMAX=0 . 

TMINMOO. 

ETAMAX=0 . 

ETAMIN=1. 

RMIN=1000. 

RMAX=0 . 

TTOT=0. 

ETAT=0 . 

LITE=0 . 

DARK=0. 

L=LENGTH 

««««#«« 

BEGIN  ITERATIONS 

DO  500  N= 1,1 TER 

SUN1=S0RT((C0S(THETA)#C0S(L0NG)-SIN(THETA)#C0S(XINC) 

+  #SIN(L0NG))«»2.+(SIN(THETA)#SIN(XINC)»C0S(ZEN)-(C0SCTHETA) 

+  #SIN(L0NG)+SIN(THETA)#C0S<XINC)#C0S(L0NG))«SIN(ZEN))»«2.1 
SUN2=ABS( C0S(THETA#PI/5 .319)) 

I F( THETA . GT . PI )SUN2=ABS( C0S(  ( 2 . #PI -THETA ) *PI /5 . 31 9 ) ) 
REST=RESTR*(  1 . +ALFA)K  TEMP-293. 15) ) 
R=4.»REST#L/CPI#D»«2. ) 

I F< cose  THETA ).GT.C0S( PI -SHADOW ))G0  TO  130 
:F((VB(N)«L/LENGTH)*i*2.  .LT.4000.#P«R)GO  TO  3000 
AMP=VB(N  )*L  ( 2 .  «LENGTH#R )-  .  5#S0RT(  ( VB(  N  )^tL/  ( L£NGTH«R ) ) 
»*2.-4000.#P/P) 

GO  TO  150 
130  AMP=-XAKP 

;F(NN.E0.2)THEN 

AMP=-(VOC-V3(N)#L  LE!iGTH-VC0'<7;  P 
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I FC  VB(  N ) itL/  LENGTH .  GT .  VOC- VCONT ) rt!>'iP=U  . 

END  IF 

I F( NN . EO . 3 ) AMP=-THRUST ' ( L*VMAB<  N ) ) 

1 50  RHO=RHOR* ( 1 . -CRHO( TEMP- 273.15)) 

MASS=RH0»PI*(D#»2. 

FORCE=AMP«VB( N ) *u'  ( LENGTH#V ) 

I F( NN . EQ . 3 . AND . COS ( THETA  > . GT . COS  ( P I - SHADOW ) ) 

+  FORCE=-THRUST 

SOLAR=SOL#A*D#L«SUNl 
EARTH=ERAD»EM«F#PI »D«L 
ALBEDO=ALB#SOL»A«F*PI  «D)tL»SUN2 
0HM=(A«P#«2.)«R 

VA I R=465 . # { RORB/ RE ) «COS ( X I NC#COS (THETA)) 
VREL=SQRT( VAIR*«2 . +V*#2 . -2 . #VA IR#V#COS( XI NC* 

+  SINCTHETA))) 

AIR=(3. 113D-1 1 )»D»L»VREL)n»3. 
RAD=PI«D»L)fB0LT«EM)KTEHP#»4. ) 

OTEMP=TEMP 

I F( COS( THETA ).GT.COS( PI -SHADOW)) GO  TO  180 
TEMP=TEMP+ ( EARTH+OHM+AIR-RAD )#DELT/ ( MASS*C ) 
DARK=DARK+DELT#AMP«VB(  N )  )tL/LENGTH 
GO  TO  200 

1 80  TEMP=TEMP+ ( EARTH+OHM+A I R+SOLAR+ALBEDO-RAD ) *DELT 

+  /(MASS»C) 

LITE=LITE+DELT»AMP#VB(N)*L/LENGTH 
200  ETA=1 .-LENGTH#ABS(AMP)»R/(VB(N)#L) 

WRITE  TO  DATAFILES 

IFCNN.EQ.DTHEN 

WRITE{51,*)T(N).TEMP 
WRITE(52,#)T(N).R 
WRITE(53.»)T{N).AMP 
WRITE(54,»)T(N). FORCE 
WRITE(55,»)T(N).ETA 

END  IF 

IFCNN.EQ.DTHEN 

WRITE(61 ,*)T(N),TEMP 
WR1TE(62.*)T(N),R 
WRITE(63.«)T(N).AMP 
WRITE( 64, #)T(N). FORCE 
WRITE<65.«)TCN) .ETA 

END  IF 

IFCNN.EQ.DTHEN 
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c 

c 

n 


c 


WPiTEc": .*;T(N) .TEMP 
WRIT£<  ?2.*)i’<N)  .P 

wp;te(  .amp 

WRITE; 74. »)T(N). FORCE 
WP!TE;^6.»)T(N) .ETA 

END  IF 

#«**#*#*»**'»*(»#****M###**»**#**HHt******* 

CALCULATE  PARAMETER  AVERAGES  i«  VARIATIONS 


IF(R.LT.RMIN)RM1N=R 
!F(P.GT.RMAX)PMAX=P 
I F ( TEMP . GT . TMAX ) TMAX=TEMP 
I F ( TEMP . LT . TM I N )TM I N=TEMP 
I F ( ETA . GT . ETAMAX ) ETAMAX= ETA 
I F ( ETA . LT . ETAM I N ) ETAM I N=ETA 
ETAT=ETAT+ETA 
TTOT=TTOT+TEMP 
THETA=THETA+RATE»DELT 


500 


c 

c 

c 

•c 

c 

c 


c 

c 

c 


L=L»( 1 .+AEXP»(TEMP-OTEMP) ) 
CONTINUE 
TAVG=TTOT/ITER 
ETAVG=ETAT/ITER 


MATCH  ENERGY  LEVELS  OF  GENERATION  &  THRUSTING 


(ITERATION  WILL  MATCH  ENERGY  TAKEN  FROM  ORBIT  WHILE 
GENERATING  POWER  IN  SHADOW  WITH  ENERGY  INJECTED  INTO 
ORBIT  WHILE  THRUSTING  IN  SUNLIGHT) 


IF(KEK.EQ.2)G0  TO  600 
IFCABScABS(LITE)-uAPK).LE.ERR)GC  TO  600 
I F ( ABS ; ^ I TE ) - DARK . LT . 0 . )THEN 
IFCNN.EO.; )THEN 

IF(XAMP.EG.PREV)GC  TO  600 
I F  (  DAP?;- ABS(  LITE ) .  GT  .  i  5 .  *ERP  )XAMP-XAFP-^  . 
IF(  DAP?;- ABS(  LITE) . GT . 5  .  # ERR.  AND .  DARK- 
ABS;  L .  TE  ; .  LE .  1 5 .  *EPP  );-;AMP=XAMP-  .  0 1 

IF;  DAP?:-ABS;  LITE) .  LE  .5 .  *ERP).\AMP=:-,AM.;-  *  .  I; 


1 


.  r ' . « .  i .  t.  .  c. ' .  h  r.  i’< 
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vv 
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IF(VOC.EQ.PREV)GO  TO  600 
I F( DARK-ABS( LI TE ) . GT . 1 5 . #ERR ) VOC=VOC+ 100. 

I F( DARK-ABSC  LITE ) . GT . 5. «ERR . AND . DARK- 
+  ABS( LITE) . LE . 15 . «ERR) VOC=VOC+ 10 . 

I F( DARK-ABSC  LITE ) . LE . 5 . *ERR ) VOC=VOCt 1 . 

END  IF 

IFCNN.E0.3)THEN 

IF(THRUST.EO.PREV)GO  TO  600 
I F( DARK-ABSC  LI TE ) . GT . 1 5 . *ERR )THRUST=THRUST+ . 1 
I FC  DARK-ABSC  LITE) . GT . 5 . »ERR . AND . DARK- 
+  ABSCLITE).LE.15.«ERR) THRUST  =THRUST+ . 0 1 

I FC  DARK-ABSC  LITE ) . LE . 5 . «ERR ) THRUST=THRUST+ .001 

END  IF 
GO  TO  90 

END  IF 

I FC  ABSC  L ITE ) -DARK , GT . 0 . )THEN 
IFCNN.EQ.DTHEN 

IF(XAMP.EO.PREV)GO  TO  600 
PREV=XAMP 

I FC  ABSC  LI TE ) -DARK . GT . 1 5 . »ERR )XAMP=XAMP- . 1 
IF(ABS<LITE)-DARK.GT.5.*ERR.AND.ABSCLITE) 

+  -DARK.LE.15.»ERR)XAMP=XAMP-.01 

I F<  ABSC  LITE )-DARK . LE . 5 . »ERR)XAMP=XAMP- . 00 1 

END  IF 

IFCNN.EQ.2)THEN 

IF(VOC.EO.PREV)GO  TO  600 
PREV=VOC 

I FC  ABSC  LITE ) -DARK . GT . 1 5 . #ERR ) VOC=VOC- 1 00 . 

I FC  ABSC  LITE ) -DARK . GT . 5 . »ERR . AND . ABSC  LITE ) 

+  -DARK.LE.15.»ERR)V0C=V0C-10. 

I F  C  ABSC  L I TE ) - DARK . LE . 5 . #ERR ) VOC= VOC- 1 . 

END  IF 

IFCNN.EQ.3)THEN 

IFCTHRUST.EQ.PREV)GO  TO  600 
PREV=THRUST 

I FC  ABSC  LITE) -DARK . GT . 1 5 . #ERR )THRUST=THRUST- . 1 
IFCABSCLITE)-DARK.GT.5.*ERR.AND.ABSCLITE) 

+  - DARK . LE . 1 5 . »ERR )THRUST=THRUST- . 0 1 

I FC  ABSC  LITE)- DARK . LE . 5 . «ERR )THRUST=THRUST- .001 

END  IF 
GO  TO  90 

END  IF 


PRINT  RESULTS  TO  SCREEN 
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60U  WRITE(6,*)  A  '  EH  -  AV  I  ■  MX  T  /  iiN  T  AV  EF  MX  £F 
-  EF 

WRI TE;  6 . 2 )  A .  EX .  TAv'G .  TMAX .  THIN .  ETAVG .  ETAMAa  .  ETAM I  fJ 

2  FORKATc  I/. ,  F3. 2 ,2>. . F3. 2 ,3FT .  i  . 3F7 . 4 ) 

WRIT£(6.o)  MAa  resistance  =  .RMAX,  OHMS 

WRITEC6.3)  MIN  RESISTANCE  =  .RMIN.  OHMS 

3  FORMAT(A.F8.3.A3 

WRITE<6.3)  FINAL  TETHER  DIAMETER  =  '.10a0.*D.  MM 
IFCNN.EQ.: )WPITE(6,3)  FINAL  CURRENT  =  .aAMP.  AMPS 

Ir(NN.EC.2)WRITE(6.3J  FINAL  VOLTAGE  =  .VOC.  VOLTS 

:F<NN.EC. j)WPITE(6,3)  FINAL  rORCE  =  .THRUST.  NEWTONS 
WRITE(6.#)  ENE.=>GV  DIFFERENCE  =  .ABSt  L;  IE)-DARK 
GO  TO  4000 

3000  WRITE(6,*5  D  IS  TOO  SMALL  TO  MAINTAIN  CON.STANT  P 
D=D+.0001 
GO  TO  90 
4000  END 


MN 
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APPENDIX  B 


Computer  Calculation  of  Variation  of  Parameters 

The  variations  in  the  parameters  of  the  orbit  of  the  tether  system  are 
evaluated  as  discussed  in  Section  3.3.  The  computer  program  which 
implements  the  mathematical  model  described  there  is  included  in  this 
appendix.  It  employs  a  fourth  order  Runge-Kutta  numerical  integration 
code  to  compute  the  changes  in  the  various  orbital  elements  with  a  time 
step  of  60  seconds. 

All  computations  are  referred  to  the  e,  p,  h  coordinate  frame, 
although  inertial  coordinates  (used  to  specify  the  terminator)  in  an 
earth-centered  reference  frame  are  also  calculated  at  each  iteration. 

The  equations  used  are  as  follow: 


^  -  r  sin  6 
dt  “  h  sin  i 

dt  '  /*ae 

di  =  r  cos  9 
dt  h 


^d 

ih  •  ^d 

[(l  •  Y][r  .  y  +  [pa 

ih  •  ^d 


^  -  zH  \  fh  (cos  f  +  e)  h)  T 
dt  =  ll  p 

The  sixth  classical  orbital  element, 
passage. 


-  r 


(p  +  r)  sin  f  /]  a^  "  cos  i  ^ 
of  course,  is  time  from  perigee 


The  components  of  the  disturbing  acceleration  vector,  with  the 
exception  of  the  force  due  to  the  earth's  oblateness,  are  as  described  in 
Section  3.3.1  of  the  main  body.  The  total  acceleration  due  to  gravity  is 
determined  as  the  gradient  of  the  disturbing  function  associated  with  the 
earth's  gravitational  field  as 


The  disturbing  acceleration  due  to  oblateness  is  that  portion  of  the  total 
acceleration  not  associated  with  normal  two-body  motion,  i.e.,  the 
infinite  sum.  Taking  only  the  first  three  harmonics,  as  described  in 
Section  3.3.1,  the  vector  components  of  the  disturbing  acceleration  due  to 
oblateness,  referred  to  e,  p,  h  coordinates,  are 


JL. 

g.e  =  ^2 

'  4 
k=2 

.r  , 

K  1 

[P|^^l(C0S  ^)cos  f  ■ 

-  P|^(cos  ^)sin  i  sin  wj 

4 

r  1 

1  ^ 

JL. 

g.P  =  ^2 

k=2 

e 

[r  . 

K  1 

[P^^j(cos  ^)sin  f  - 

-  P|^(cos  ^)sin  i  cos  wj 

4 

K 

1  \ 

JL. 

g.e  =  ^2 

k=2 

e 

.r  , 

P|^  (cos  ^)cos  ij 

These  are  the  components  used  in  the  computation. 


ooooooo  onoooof) 


C  THIS  PROGRAM  COMPUTES  THE  VARIATION  OF 

C  INITIALLY  SPECIFIED  ORBITAL  ELEMENTS  FOP  A 

C  TETHER  SYSTEM  OPERATING  IN  A  SPECIFIED  MODE 

C 

C  DIMENSION  ARRAYS 

f' 

c 

COMMON  I ,W(20000). INC(2000C) .LON6(200UG).R<2000U.3) . 

+  VC  20000,3) .FC 20000, 3) ,TF( 20000 ) .A( 20000 ) .H( 20000 ) .P,RD,>. . 
*  Y,2,E(1),EC(20000).T(20000),DELT.22.D(3),Q1,02.03 
REAL  INC, MU, LONG 
INTEGER  I 

DEFINE  SYSTEM  PARAMETERS 

#«*»#«*«#»»#«#««««««««*«*»*««*«««««*««*« 

PARAMETERC  HEI GHT=400 . D3 , RE=6378 . D3 , MU=3 . 986D 1 4 ) 

PARAMETER ( P I =3 . 1 4 1 5927 ) 

«#««««««#«#««««««#««««««#»«««#««««»»#«#« 

SPECIFY  INITIAL  CONDITIONS 


11  =  1 

DELT=60 . 

ZZ=-1. 

A(1  )=RE-^HEIGHT 
T( 1 )=0. 

INC( !  )=.4<>74 
LONGC 1 )=0. 
WC1)=PI  2. 

El '  1  3  =  .  o  i 


n? 


f'/oon  ooooooo  ooooooo  ooooooo 


OPEN  DATAFILES 

*4***#**#*  *»***#*««'* 

0PEN<21 .FIL£=  INC  .STATUS=NEW  ) 

0?EN(22.FiLE=  ASC  .STATUS=  NEW  ) 

0P£N(23.FILE=  AXIS  .STATUS=  NEW  ; 

0PEN(24.FIL£=  ECC  .STATUS=NEW  i 
0PEN<25.FILE=  LONG  ,STATUS=  NEW  5 
0PEN<26.FILE=EXP  .STATUS=  NEW  ) 

NUMERICALLY  INTEGRATE  VARIATIONAL  EQUATIONS 

*«»««*«4(«*»»»#**«««««**««*#*4«**#**#***4 

DO  1=1,12000 

»*««««««««•«««*«««««««**«***»«««*«*««««* 

COMPUTE  INERTIAL  COORDINATES  OF  TETHER 

»«««*«««»««*««»«««#«#«#*«««««««*#»«»*««* 

P=(H<I)##2.)/MU 

RD=P/ < 1 . +EC( I ) «C0S( TF( i ) ) ) 

IFCRD.LE.REIGC  TO  500 

X=C0S(TF( I )+W( I ) )*C0S( LONGC I ) )-SiN<Tr ( i  )+WC i ) ) 
+  »COS(INC(I))#SIN(LONG(n) 

Y=SIN(TF<I)+W(  1 ))*COS( INC(  I ) )#C0S( LGNG( . );♦ 

+  COS<TF(I)+W(n)#SIN(LONG(I)) 

2=S:N(TF( I )+W( I ) )*SIN( INCC I ) ) 


CALCULATE  ORBITAL  ELEMENTS 

lALL.’  DEFINE  STATE  (POSITION  A  VELOCITT)  /tCT  .r 


II3 


n  Art  icrt  A.t  fwi 


•SSi 

m 

I 


>«*•«»*»****•* 


CALL  STALE; Lf! I ) ,EC(  I ) ; 


*»***»»»*»»*»*#(***###»»##)►*##»*#**#*##♦* 


CALCULATE  DISTURBING  ACCELERATION  AT  TIME  INSTANT 


#<►»*##»»*♦♦•*###«###«»*#***•»#*<►«»#**#*#* 


CALL  ACC(TF( I ; .EC< I ) ) 

T( 1*1  )=T( I )*DELT 

»**»#####*#«#«###»«»###»##«#####*♦•*«#♦»• 

LONGITUDE  OF  LINE  OF  NODES 

LONG;  I  ♦  I  J=L0NG(  I  /  *?D*DELT’‘SIN(W(  I  )-'rFc  I  )  >*F(  1 ,3 
•SINC iNC( I ) ) ) 


sec;  ■'■a  :o?  axis 


ANGULAR  MOMENTUM 

H< I+l )=H( I ;+DELT#(R( I , 1 )»F( I ,2)-R( I ,2)*F( I , 1) ) 

TRUE  ANOMALY 

»##»**###»«#»»#*<(#*«##«#»**»##*#»♦###### 

CALL  TRUE(DELT,EC(I).TF(I+1)) 

ORBITAL  INCLINATION 

INC( I+l )=INC( n+DELT#RD»COS(W( I )+TF( I ) )*Fi I ,3) 

+  /H(I) 

ARGUMENT  OF  PERIGEE 

W< I+l )=W( I )-(P#COS(TF< I))«(R( I ,1)#F( I , 1 ;+R( I .2) 

■*  (F(  ;  ,2);'(P+RD)»SIN(TF(I))«(R(I.1)«F(I.2)-RCI.2)«F(I.1))) 
*  +DELT/ (RD+ECC I )*H< I ))-COS( INC( I ) )»(LONG( I+l )-LONGC I ) ) 

.1.:^-'  ,Nl  CuPPE't.  L-Pr-v^iiCN 


,  F (  V  !*  ♦  F *^2 .  ) .  : 


. .  .  U  .  ^  .  '1M'< 


cc?^?:jTE  vA?:Ar:r;:;  in  cp3;:a:  altiTuLIS 

r iGEr.  —  ^A'  1  /  >*(  1  .“ lC .00*^. 

A?0GEE=(A( I )*( I .+EC( I ) )-RE)  lOUU . 

GO  TO  600 

WR I TE(6.*)1  .-RADIUS  BELOW  SURFACE  OF  EAR,r: 
bC  i G  6000 

♦  •in*##**##**-#*#####*###***-*#*##**  ♦«#*♦#■»• 

WRITE  TO  DATAFILES 

»»»##»*##»#»**■###»#»*#*«*#»♦»#***♦*<►**** 

IF<n.£0.5)THEN 

WRITE(21 .»)T(I)/86400. .INC(I)#180./PI 
WRiTE(22,«)T( I )/86400. ,W( I )*180./Pi 
WRITE(23.»)T(  I  V86400.  ,A(  I  )-  1000. 
WPiTEc24.#)T( I )/86400. ,EC( i ; 

WRITER  25, *)T< I )/ 86400 . .LONGC I )*180 ?: 
WRiTE(  26,  <♦)!(  I  )  86400  .  , PER i  GEE 
GO  TO  5000 


GO  ;C  bOOO 


lND  if 


AT  THr  INo“AN:ANEnL;S  P^IM 


m 


•ffi 


1 


Cur^MGN  .  .W:2'J000  5 . 1NC( 20000  ). L,0!iG< 200CG  ;  .Pc  20000  ;  , 

V: .  :;0u:. .  : ; .  r  c  20000.3)  ,TF(  20000 )  .Ac 20000  ), He  20000  )  .P.PD.X. 

'.  . ..  .KC :  ;  .  ECC  20000  )  ,TC 20000 )  ,DEi.T.Z2  ,  DC  3)  ,  Cl  .02 .03 

'EAL  : :;c . '12 , lomg . l i b . lat . mass . :. . l i bmag 

2‘<7LGE.^'  1 

::MEfis:GN  sc 3; .currcs) .vpelcS) .A:i?(2oooc..?) .sclarcsoooo.s; 

2 1  MENS  I  ON  :.IBC2COOC,0-i;,EMF(2000G.3>,GRAv(20000.3; 

PARAMETER C  L=2 . D4 . AMP=3 . . TEMP= 1100., MASS=25 , D3 . R£=63T8 . D3 ) 
PARAMETER c  COEFF= 1 .25, FLUX=4 . 4D-6. Pi=3 .141 5R27 , DkAG=2 . 2 ) 
PARAMETERc  Mu=3 . 986D1 4 . DI AM= . 003 ) 


COMPUTE  MAGNETIC  FIELD  STRENGTH 


*#»»#»#»#»##»»####)♦««»##*#»###<*#######*# 


«*Iy 

**tS' 


n 


CALL  MAG(TRF) 

BQ=D< 1 )*COS<  LONGC I ) )+D{ 2)«SIN( LONGC 1 ) ) 

BR=-D( 1 )#SIN( LONGC I ) )+D(2)#C0SC  LONGC I ) ) 

Be  1 )=BO»COSCWC I ) )  +  (BR»COS( INCC I ) )  +  Dc3)«SlNC INCe I ) ) l*SINCWc I ) ) 
BC2)=-BQ»SINCWC I ) )  +  CBR«COS( INCC I ) )  +  DC  3)#SINC INCC I ) ) ) 
*COSCW<I)) 

BC3)=-BR«SINC INCC I ) )+DC3)«C0SC INCC i )) 


)♦*##***# 


COMPUTE  CURRENT  VECTOR 


.  •  n  .  .. 


t,Mr;  .  ..  '  .... r'c  ;  : 

EMF<  !  ,  '  J  =  ZF  * L* (  B(  ii  V  .  ■  '.  ;  ♦..Jr'F  >  ^  -  - 


.  ERMAG=-»6do  .  *  V..  V,-  3( LA i ) #PD-  PE 
i  Ft  ]  .Em  )7HE!, 

VREL; ; ;^-VERMAG*( t -SINtPSI )#COSt  LONGC I  )  >GG3t  PG: ; 

•  *siNt LONGt  1 ; ;  )*cos(W(  I )  )-‘(ccs(Psi  )#cost  LCiiGt  I )  ;-s;Ntps;  ■ 

-  *s:n( LONGt ; ; ; )#cos< iNCt ; ) )«sin(W( d)) 

VRELC2)=SQRT(MU.  At  1 ) )-VERMAG*t (SINCPSi )*CUb( LONo 

*  ( I ) )-C0S(PSi )*SIN( LONGt i ) ))*SINtWt I ) )+tCCScPSI )*CCS 

-  t LONGt I ) )-S!N(PSl )#SINt LONGt I)))#COStWt I ) )#COSt iNCt i ) ) t 
ELSE 

VRELt 1 )=Vt I ,l)-VERMAG«tt-SINtPSl )*CDSt LONGt i 
^  COStPSI)*SINtLONGtl)))»COStWtl))+tCOStPSl )#COStLONGt I ))- 
+  SINtPSn#SINtLONGtI)))»COStINCtI))#SINtWt  I))) 

VRELt2)=Vt 1 .2)-VERMAG#t  t SINtPSI )*COSt  LONGt I ) )- 
+  COStPSI)»SINt  LONGt  I ) )  )itSINtWt  I )  )  +  tCOStPSI  )»COSt  LONGt  I )  )- 
+  SINt  PSD  *SINt  LONGt  I ) )  >«COStWt  I )  >»COSt  INCt  I) ) ) 

END  IF 

VRELt  3)=-VERMAG«t  SINt  PSI )»SJNt  LONGt I ) )-COSt  PSI J  *COSt  LONGt I ) ) 

-  #SINtINCti)) 

REL=SQRT t  VRELt  1  )»«2.  +  VRELt  2)**2 .  ■►VRELt  G )  **2  .  ) 

COMPUTE  AIR  DENSITY 


C  'J  CO 


“..'■'  .  .  -  .  ■:)*Dr"’AG>iJ  1  A!'!*l*COSv  "i-t  A  ) »DENh)'*.’^£^* ''RE-- '  J  ' 


»*»*-i»»**)^*i* 


COMPUTE  TETHER  NORMAL  AREA 


?ERr=SQRT( cCOS(ALFA)»COS(W(i )tTRF)«CGScLONG( I ) )-SIN( ALFA)* 
;:.'vWC  I  )  +  TRF5»COS(LONG(I  ))-COSCALFA))*SIN(W(I)+TRF)» 

OS( INC( I ) ;«S:N(L0NG( I))-SIN(ALFA)*COS(W( I )+TRF)*COS( iNC( I))* 
*  SiNCLONGC  I ;  n»*2.-KC0S(ALFA)«SIN(W(  I  )+TRF )«SIN(  INCt  I  ))  + 

-  S:n<ALFA)#CCS(W(I)+TRF)«SIN(INC(I)))*«2. ) 

•  *«*««*«*  •)* -im  It 

CALCULATE  SOLAR  RADIATION  ACCELERATION  VECTOR 


■  fRL*S0RT(X»*2.*Z**2. ) .LE.RE.AND.Y.LT.O . )THEN 

Ju  vO  ~  i  ,  O 

S0LAR( I .J)=0. 

END  DO 
-br, 

: . :  ;  =  -tccEFF#DiAM*:.-?E:p?^siNc  :NCf :  - 
:i:c  ;  n.cob(W( : ; ;  ■ 


.  j  .  .  «  \  i  .  I V .  ■ 
.'.Tv  .  '  1  ...  . 


L I BMAG— 3 (2.*? / 

LIB<  I  ,n=-LIBMAG«<COS(TRF;-^ALFA«SlN(TRi-  ;  ) 

LIB( I .2)=-LIBMAG»(SIN(TRF;-ALFA#C0S(TPF;  ; 

LIB( 1 .3)=0  . 

*»»«#«#*»•»♦#»»#»»»*##««*##»**#***♦*♦**»* 

CALCULATE  TOTAL  DISTURBING  ACCELERATION  VECT 

DC  j=; .3 

:< :  .U)=AIP( : .:;-GPAV; ; 


.j)^l:b 


.J)-SGLAP 


oofiooooonri 


SUBROUTINE  TRUE( DT.ECCC.TRF) 


COMMON  I . W( 2000U ) , INC( 20000 ) . LONOC  20000  ) . Rl 2O000 , 3 ) , 

+  V( 20000, 3). F( 20000, 3). TF(20000).A( 20000 ).H( 20000 ),P,RD.X, 
+  Y,:,E(  1 )  ,EC(  20000  ),T(  20000  ).DELT,:Z,Ii(  3)  ,Ql  ,02,03 
REAL  INC, MU, LONG 
INTEGER  I 

PARAMETERC  MU=3 . 086Dl  4 ,  PI  =3 . 1 41 5<527 ) 

TRF=TF(  I  )+DT*(h<  i  )/RD*^*2.'t(P«C0ScTF(  i  );*(R(  I ,  i  ;*F(  1 , 1)  + 

+  R< I ,2)mF< I ,2)-(P+RD)»SIN(TF( i ) )*(P( 1 . 1)»F(  1 ,2)-R(  I  ,2)« 

^  F(  I  ,  1  )))^(RD*ECCC<fH<  i ) ) ) 

RETURN 

END 


COMPUTE  TETHER  POSITION  &  VELOCITY  VECTORS 

SUBROUTINE  STATE (TRF.ECCC) 


COMMON  I.W(20000),INC<20000),LONG(20000),R(20000,3), 

+  V(20000,3),F(20000,3),TF(20000),A(20000),H(20000),P,RD,X, 
+  Y,Z,E(1),EC(20000),T(20000),DELT,22,D(3),01 ,02,03 
REAL  INC, MU, LONG 
INTEGER  I 

PARAMETER( MU=3 . 986D 1 4 ) 

R( I ,l)=RD#COS(TRF) 

R(I ,2;=RD*SIN(TRF) 

R(I,3)=0. 

VR=MU»ECCC»SIN(TPF)/H( I ) 

VT=MU*(  1  .+ECCC«CC.SC7PF;  )/H(  I  ) 

V( : . 1 ;=VP#COS<TPF;- VT^SINCTRF) 

V  C  .  ,  w;! ' "  v'P-^G i  N f  i  r’:  ?  *  .'  i  C  THi'  / 


>  _ 


This  SSb?.-Si  :T,t  CAiaiArES  :HE  F:ELD  STPr.F.u:-'  ^  £Ar  ;-  . 
MAGNETIC  -lELT  (MODELLED  AL  AN  UCTOFOLE)  AT  Tnr  '.CLAN-a'.-  , 
POINT  CN  THE  .o;'.LA:iNG  CFBIT 


SUBROUTINE  XAGv  TRE 


'^■oMMU.  1  .  ,  *  '  ^v. # , .  .4*^.  ^  V. - .  ^  ,  4^*,iNG  ^ <  '  .r'  i_ U . 

V  (  **^0 0  ■*  ^  (  ’ ■  ‘  ^  f"'  rT  ^  ^  JT  '  '  *  , ,  '  '  ' ' 

I  ,i.,t(  1  /  .^'''i^UGOoi.  i(DCGU-  ,  '■  i  .01  ,  'JL  , 

^  i  .M E.N^  i  i..N  r”r '  t ,  *4 PO i  -t » ^  T ,  *1 '  .  o ( *4 ,  n  '  . . ' i  4  ,  *4  ,  rM  ■  4,4- 
DIMENSION  BB(  3)  .GG(4.4;  .hr.:  4.4) 

REAL  INC .  .Mu  .  LONG .  M  . N  .  -.C .  MLo.NG 
INTEGER  I 

PARAMETEPc  RE=6378 .  D3 . RGTPATE=  c .  3]  2D-5 .  P I  =3 .  I  4 1 5^.,  “ .  UPDEP=4  . 
DC  K=i .ORDER 

DC  J=l. ORDER 

G<K,J)=0. 

HM(K.J)=0. 

END  DO 

END  DO 

G(2,n=-30186.D-A> 

G(2.2)=-2036.D-9 
G<3,l)=-1898.D-9 
G(3,2)=2997.D-9 
G<3.3)=1551 .D-9 
G(4,1)=1299.D-9 
G<4,2;=-2144.D-9 
G(4.3)  =  l296.D-<3 
G(4,4)=805.D-'^ 

HM(2,2)=5735.D-‘5 
HM(3,2)  =  -2124.D-«J 
HM(3.3;  =  -3''.D-9 
HM(  4,L;=^-361  .L-^ 


3(N, ;  ;=S(N-i . ,  1  )<*c  c 2.'*(N- 1 .  )-i . ) /  (.’■i-i . ) ; 

DC  ORDER 

IF(M.EG.2. )THEN 
DEL-1 . 

C.LSE 

DEL-C . 

END  ;f 

.  )#SQRT(ABS(  (N-M+1 .  )«(DEL+i  .  ) 


END  DC 
DO  J=1  . 3 

3B  V  ^ / —0  . 

END  DC 
ANG=TRF+W( I ) 

THETA=ACOS( SIN< ANG)*SIN( INC( I ) ) ) 

IF(ANG.EG.PI/2. .0R.ANG.E0.3.«PI/2. )THEN 
ALFA=ANG 

ELSE 

ALFA=ATAN(COS( INC( I ) )«TAN(ANG) ) 

END  IF 

IFCANG. GT. PI/2.. AND.ANG.lt. 3. »PI/2.)ALFA=ALFA+PI 
iF(ANG.GT.3.<*PI/2..AND.AN6.LT.2.#PI)ALFA=ALFA+2.»PI 
MLONG=ALFA+LONG< I )-T( I )»ROTRATE 
DO  N=2, ORDER 

DO  f1=l  .ORDER 

IF(N.EO.M)THEN 

?S(N,M)=SIN<THETA)*PS(N-1 . ,N-1 . ) 

ELSE 

PS( N,M)=COS( THETA ) *PS c N . .M)- 

♦  KC(N.M)»PS(N-2. .M) 

END  IF 


END  DC 


• 


I 


’mm. 


END  DO 

END  DC 

DO  N=2, ORDER 

DO  M=l. ORDER 

GG(N,M)=S(N.M)«G(N,M) 

HH(N,M)=S(N.M)*HM(N.N) 

END  DO 

END  DO 

DO  N=2, ORDER 
SUM  1=0. 

SUK2=0 . 

SUM3=0 . 

DO  M=1.N 

SUM1=SUM1+(GG<N.M)»C0S((M-1 . )*MLONG)+HH(N.M)# 
T  SIN((M-l.)*MLONG))»PS<N,M) 

SUM2=SUM2+(GG(N,M)*C0S((M-1 . )*MLONG)+HH(N,M;« 
+  SIN<(M-l.)#MLONG))*PP(N,M) 

SUM3=SUM3+(M-1 . )#(-GG(N,M)»SIN( (M-1 . )*ML0NGn 
+  HH(N,M)«C0S(<M-1 .)»MLONG)#PS(N,M) 

END  DO 

BBC  1  )=BB(  1 )+( (RE/RD)<HKN+1 . )  )«N#SUM1 
BB(2)=BB(2>+< (RE/RD)##(N+1 . ) )»SUM2 
BB<3)=BB(3)+<<RE/RD)«*(N+1 . ))«SUM3 

END  DO 
Bfi(2)=-BB<2) 

BBC3)=-BB(3)/SINCTHETA) 

DC l)=(BB(l)»C0S(PI/2.-THETA)+BB(2)*SIN(PI/2. -THETA) )« 

+  C0S(ALFA)-BBC3)icSINCALFA) 

DC2)=CBBC1 )«C0SCPI/2. -THETA)+BBC2)*SINC PI/2. -THETA ))» 
SINCALFA)+BBC3>»C0SCALFA) 

DC3)=6BC 1 )»SINCPI/2.-THETA)-BBC2)*C0Scpi/2.-THETA) 

RETURN 

END 


SCBRCl  7 ; ‘.E 


SUBROL'Tl.N’E  EARTH (TRF) 


COKMUN  :  ,  W(  ) ,  iNl..'  l^JUOO  ) ,  uuhGC  iiOJiJi.-  -  ,  . 

V{200CC.B'/  .FC2UUOO.J)  .rr(200C0  )  .A(2000C  ).  H(  20000  ).  P ,  R 
-  Y,Z,E(  i  KEC(2U00G>  ,T(200Q0)  .01  .02.. 

DIMENSION  ZZZ(3)  ,:.P(5)  .LPP(6)  .GRAVGos; 

REAL  INC.MU.LONG.LP.LPP.OJ 
INTEGER  ; 

PARAMETER(  Mu=3 . 986D1 4 . RE=63T6 . D3 ) 

JJ( 1 )=C  . 

JJ(2)=1 .08263D-3 

JJ(3)=-2.54D-6 

jJ(4)=-l .61D-6 

22Z(  1 )=SIN( INC( I ) )*SIN(w( i ) ) 

Z2Z<  2)=SIN( INC( I ) )*COS(W( I ) ) 

ZZZ<3)=C0S( INC( I )) 

XI=SIN(W(I )+TRF)mSIN(INC( I)) 

LP( 1 )=1 . 

LP(2)=XI 

LP(3;  =  .5#<3.«XHn<2.-l . ) 

IFCXI .LT.O. )THEN 

LP(4)=.5»<-5.#<ABS(XI))«»3.-3.*XI) 

ELSE 


LP(4)=.5*(5.#XI««3.-3.i(XI) 
LP<5)=<35.*<ABS<XI))#»4.-30.«XI»#2.+3. >/8. 

LPP(1)=0. 

DO  XK=2,6 

LPP(XK)=XI»LPP(XK-1 . )+(XK-l . )»LP(XK-1 . ) 

END  DO 
DO  J=l,3 

SUM=0 . 

DO  XK=2,4 

SUM=SUM+JJ(XK)»< ( RE/RD) »«XK)*(LPP(XI>2. 
-  »R(  I  ,  J).-  RD-LPP(XKM  .  )*2ZZ(J)) 

END  DC 


oooooooooo 


♦  ♦  ♦ 


c 


CALCuLAiE  CHANGE  .N  ECCEN. RICH') 


SJBRCi. 


iERlvc 


COMMON  ■  ,  '  /LuOOO  ) ,  INCc  OOCCC  J ,  l,CNG(  20C0C  )  .  P\  ZJOOC  .  -r  /  . 
v<2G0UtJ  :  HCCGCi  .3>  .7F;  ^CGCC  )  .ri<  200G0)  .  nHGGGG  )  .r  .RL..-. . 
*  V  ,EH  )  .ECC2G00C)  .TC-COG;; :  .  DELT . ZZ  .  Cc  3 )  .  o  .  .02  . ■ 

REAL  INC.M..  ,_ONGi 
INTEGER  : 

PARAMETER  cm;: =3.^58621 4; 

0E=  ( ( P^A  c  .  '-.''2**2.  j  )* .'  '•  .  ,i  )+»<  .  ,2/*fC ,  ,t.'  /  ■^cRCl  ,i 
+  V(I,l)+R(  4  .2/»’»'c  I,2))*cR'  i.I  /  »FCi,l)i'R(  I  .2)*F(  i  MU* 

^  A< I )*EC( I )  ) 


RETURN 

END 


***»*«*»*»#»*»*###«»#*#»##*###*•###*#»#*»##«#***»»*#*****«#***## 
RUNGE-KUTTA  INTEGRATION 
«»##«»*#»«*##»#*»*##»«»##**##*##«##***»* 

SUBROUTINE  RK(TRF,ECC) 

COMMON  I ,W<  2GG0C ) , INC( 20000 ) . LONGC 20000 ) ,R( 2C00C . 3) . 

*  V( 20000, ) .FC 20000, 3) ,TFC 200CC A( 20000 ). H( 200C0 ) ,P.RD,X. 

1 1 L.  ,c.C  1  /  ,E*4-C4i.0./v-'.’' ,  I  ( i^ODOO  ) ,  t-'L .  .  .4...  ,DC3)  .ui 
REAL  I  liC ,  LIjNG  . ;  ,  i,2 .  K3 .  K4 
INTEGER  I 


12(. 


•  -  o' 


APPENDIX  C 


Development  of  the  Variation  of  Parameters 

To  evaluate  the  effects  of  an  external  force  on  the  orbit  of  a 
satellite  it  is  useful  to  write  Newton's  second  law  as 

d^r  ^  n  ^ 

^  %3  E  =  j53  8r 

where  aj  is  the  vector  representation  of  the  total  disturbing  acceleration 
due  to  all  forces  except  the  primary  acceleration  associated  with  the 
gravitational  field  of  the  main  body,  and  Rj  is  the  scalar  disturbing 
function  associated  with  the  disturbing  body  mj. 

Lagrange  fully  developed  the  method  of  the  variation  of  parameters  to 
study  the  problem  of  perturbations  of  comets  moving  in  elliptical  orbits. 
Expressing  the  disturbing  acceleration  as  the  gradient  of  the  disturbing 
function,  the  result  of  his  analysis  is  the  equation 

L  ^  =  |2?P  (1) 

-  dt  laol 

where  a  is  a  six  dimensional  vector  composed  of  the  orbital  elements,  and 
L,  the  Lagrange  matrix,  is  six-dimensional  and  skew-symmetric  and  equal  to 


|9r| 

9v 

Ida 

9a 

Idal 

If  a  is  expressed  in  terms  of  the  classical  orbital  elements  as 

2^  =  in  1  tf  a  e  XI 

where  0,  i,  and  w  are  the  three  Euler  angles  as  shown  in  Figure  13,  a  is 
the  semi-major  axis,  e  is  the  orbital  eccentricity  and  X  =  -nr  the  time  of 
pericenter  passage,  then  the  elements  of  L  may  be  determined  as  detailed 
in  Reference  32,  and  the  derivatives  of  the  orbital  elements  solved  for  to 
produce  Lagrange's  planetary  equations,  as  follow: 
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nab  sin  i 


t  nab  sin  i 

^  .  -cos  i  ^ 
dt  nab  sin  i  0i 

^  =  2_  ^ 

dt  na  5X 


nab  sin  i  Qu 
b  dR 


-b 

b2 

aR 

na^e 

du 

4 

na  e 

ax 

-2 

aR 

b^ 

aR 

na 

aa 

4 

na  e 

ae 

Obviously,  L  is  non-singular  if  e  is  neither  zero  nor  one  (since  b^  = 
la^(l-e^)l),  and  i  is  not  zero. 

Gauss  eliminated  the  restriction  of  the  representation  of  the 
disturbing  acceleration  as  the  gradient  of  the  disturbing  function  by 
expressing  Equation  1  as 

da  |9r|T  ,  , 

t  df  ■  lefl  ^ 

Gauss  then  computed  the  matrix  coefficients  of  aj  to  obtain  the 
appropriate  variational  equations.  The  results  of  this  analysis,  again 
detailed  in  Reference  32,  are  Gauss'  planetary  equations 

^  _  r  sin  0 
dt  ’  h  sin  i  °dh 

^  _  r  sin  8 

dt  ■  h  ^dh 

K  =  4h  [-  p  f  pdr  *  <p  *  p  ^dt]  -  risrf  pdh 

5t  '  T  [p  f  +  fi 

^  =  -i-  [p  sin  f  a^^  +((p  +  r)  cos  t  ♦  rejaj 
S  cos  f  -  2  re)aj,  -  (p  +  r)  sin  f  a^J 


Obviously,  there  are  again  difficulties  if  the  orbital  eccentricity  is 
small  or  if  the  inclination  of  the  orbit  is  near  zero.  Of  course,  this  is 
not  a  surprise  as  the  line  of  nodes  does  not  exist  for  an  orbit  of  zero 
inclination,  and  the  line  of  apsides  is  meaningless  for  an  orbit  of  zero 
eccentricity.  A  means  to  avoid  these  singularities  occurs  if  orbital 
elements  are  defined  from  combinations  of  the  six  classical  ones  listed 
above  which  do  not  depend  on  either  the  line  of  nodes  or  the  line  of 
apsides.  For  this  analysis,  we  will  stick  with  the  classical  elements  and 
artificially  avoid  the  singularities  which  occur  when  a  particular  element 
becomes  undefined. 

Vectorial  forms  of  the  variational  equations  may  be  obtained  by 
rewriting  Equation  2  as 

da  T 

df  "  -  1^1  -d 

where  P  is  called  the  Poisson  matrix  defined  as 

P  =  -  L-1 


Following  substitution,  this  yields 


do  53 
dt  "  8v  -d 


which  may  be  used  to  directly  produce  variational  equations  of  vector 
orbital  elements  which  are  independent  of  the  coordinate  system  in  which 
the  components  of  aj  are  expressed.  With  0  composed  of  the  classical 
orbital  elements  as  described  above,  the  variational  equations  for  these 
orbital  elements  become 


da 

2a^ 

dt 

>1 

( ^ 

II 

‘  ^d 

r  sin 

■T  Ih  •  3d 

dt 

h  sin 

r  cos 

-  Ih  •  3d 

dt 

h 

dtf 

df 

^  -  COS  i 

dt 

■  ■  3v 

dt 
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h  ^  9f  , 
Ht  2  av  -d 


de  1  r 

Jt  =  I  12(!!  •  5d>  ^  -<i  •  -  (r  •  5<i. 

where  the  perturbative  derivative  of  the  true  anomaly  is 

rT-  (p..)  sinf/] 

and  the  scalar  derivative  of  the  eccentricity  is 

S  [(r  •  »)(r  •  5d>  ^  -  '■^Xi  •  5d>. 

These  are  the  equations  which  may  be  used  to  evaluate  the  changes  in  the 
satellite's  orbit  caused  by  the  presence  of  accelerations  other  than  the 
primary  one. 


^  ^  l 
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